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Abstract 

A key problem in making precise perturbative QCD predictions is to set the proper renormaliza- 

■ tion scale of the running coupling. The conventional scale-setting procedure assigns an arbitrary 
range and an arbitrary systematic error to fixed-order pQCD predictions. In fact, this ad hoc 
procedure gives results which depend on the choice of the renormalization scheme, and it is in 

r"! ■ conflict with the standard scale-setting procedure used in QED. Predictions for physical results 

should be independent of the choice of scheme or other theoretical conventions. We review current 
ideas and points of view on how to deal with the renormalization scale ambiguity and show how 
^ ' to obtain renormalization scheme- and scale- independent estimates. We begin by introducing the 

renormalization group (RG) equation and an extended version, which expresses the invariance of 
1/^ I physical observables under both the renormalization scheme and scale-parameter transformations. 

■ The RG equation provides a convenient way for estimating the scheme- and scale- dependence of a 
CN . physical process. We then discuss self-consistency requirements of the RG equations, such as reflex- 

ivity, symmetry, and transitivity, which must be satisfied by a scale-setting method. Four typical 
scale setting methods suggested in the literature, i.e., the Fastest Apparent Convergence (FAC) 
criterion, the Principle of Minimum Sensitivity (PMS), the Brodsky-Lepage-Mackenzie method 
(BLM), and the Principle of Maximum Conformality (PMC), are introduced. Basic properties 
and their applications are discussed. We pay particular attention to the PMC, which satisfies all 
a ■ of the requirements of RG invariance. Using the PMC, all non-conformal terms associated with 

the /3-function in the perturbative series are summed into the running coupling, and one obtains 
a unique, scale-fixed, scheme-independent prediction at any finite order. The PMC provides the 
principle underlying the BLM method, since it gives the general rule for extending BLM up to 
any perturbative order; in fact, they are equivalent to each other through the PMC - BLM cor- 
respondence principle. Thus, all the features previously observed in the BLM literature are also 
adaptable to the PMC. The PMC scales and the resulting finite-order PMC predictions are to high 
accuracy independent of the choice of initial renormalization scale, and thus consistent with RG 
invariance. The PMC is also consistent with the renormalization scale-setting procedure for QED 
in the zero-color limit. The use of the PMC thus eliminates a serious systematic scale error in 
perturbative QCD predictions, greatly improving the precision of empirical tests of the Standard 
Model and their sensitivity to new physics. 

PACS numbers 12.38.Bx, ll.15.Bt, ll.lO.GH 

Keywords Renormalization Group, Renormalization Scale, BLM/PMC, QCD 
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1 Introduction 

Quantum chromodynamics (QCD) is believed to be the field theory of hadronic strong interactions. 
Due to its asymptotic freedom property [1, 2], the QCD running coupling becomes numerically small at 
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short distances, allowing perturbative calculations of cross sections for high momentum transfer physical 
processes. In the perturbative QCD (pQCD) framework, a physical quantity (p) is expanded to n-th 
order in the QCD coupling as{fir)] i-e.. 



n 



p„ = Co a^M + E Ciif^r) a^s^'iftr), {p > 0) 



(1) 



i=l 



where Cq is the tree- level term, Ci the one- loop correction, C2 the two-loop correction, etc., and p is 
the power of the coupling associated with the tree-level term. The renormalization scale niust be 
specified in order to obtain a definite prediction. The calculation of the coefficients Ci{fir) involves 
ultraviolet divergences which must be regulated and removed by a renormalization procedure. The 
infinite series Pn^oo is in principle renormalization scheme and renormalization scale independent be- 
cause of renormalization group (RG) invariance [3, 4, 5, 6, 7]; i.e., the physical predictions of a theory, 
calculated up to all orders, are formally independent of the choice of renormalization scale and renor- 
malization scheme. However, at any finite order, the scale/scheme dependence from as{pr) and Ci{fir) 
do not exactly cancel, leading to renormalization-scheme and renormalization-scale ambiguities. Such 
ambiguities are well-known [8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20]. A guiding principle for 
resolving such problems is that physical results must be independent of theoretical conventions. 

It should be recalled that there is no ambiguity in setting the renormalization scale in quantum 
electrodynamics (QED) at any finite order. Mass renormalization is straightforward in QED. Due to the 
Ward-Takahashi identity [21], the divergences in the vertex and fermion wavefunction corrections exactly 
cancel, and the remaining ultraviolet divergence Z3 associated with the vacuum polarization insertions 
defines a natural scale for the running QED coupling aem(Q'^)- For example, the renormalization scale for 
the electron-muon elastic scattering due to the one-photon exchange skeleton graph in the conventional 
Gell Mann-Low (GM-L) scheme [22] is simply equal to the momentum transfer squared t = carried by 
the photon propagator. The running coupling a{t) of the Gell Mann-Low scheme is defined empirically 
from the Coulomb scattering amplitude ZiZ2a{t)/t for infinitely massive test charges. 

The renormalization scale /x^ = t = is independent of the choice of initial renormalization scale 
to since in QED 



which sums all vacuum polarization contributions, both proper and improper, to the dressed photon 
propagator. Equation (2) shows explicitly that when summed up to all orders, the result is independent 
of any initial renormalization scale to, in agreement with the RG invariance. 

The renormalization scale in QED can be determined unambiguously in any scheme, including 
dimensional regularization; the scale for different schemes are connected to the GM-L scale by com- 
mensurate scale relations (CSRs) [23], a topic which we will discuss below. The resulting perturbative 
prediction is then scheme- independent. The computation of higher-order {/3j}-f unctions for the RG 
equation is important for perturbative calculations at high orders [24, 25, 26, 27]. 

The scale-setting question is much more complicated in QCD due to its non-Abelian nature. Unlike 
QED, where there is a preferred (GM-L) scheme and a precisely known value of the coupling at zero 
momentum scale (the fine-structure constant aem(O) — 1/137.036 ■ ■ ■ [28]); in the case of QCD we do not 
have a uniquely preferred scheme and well-determined value for the coupling in the perturbative region. 
Consequently, in QCD, the uncertainty from the choice of renormalization scheme and scale must be 
treated with great care. It should be noted, however, that pQCD reduces to Abelian theory in the zero- 
color Nc — limit [29]. This analytic limit provides an important constraint on the renormalization 
scale problem in QCD. 





(2) 



i-n(t,to) ' 



where 
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In the standard procedure for a first estimate of the physical observable, one chooses a renor- 
mahzation scheme with an initial renormalization scale fir = A^™^* in E<1-(1); and then applies some 
scale-setting method to improve the pQCD estimate. After scale setting, the perturbative series for the 
physical observable (1) can be rewritten as 

n 

p„ = Co a^M) + T.C^iK) cC\K). (P > 0) (3) 

i=l 

where the new leading-order (LO) and higher-order scales and /i* are functions of the initial renormal- 
ization /i™'*, depending on the choice of the scale-setting method. At the same time, the new coefficients 
are changed accordingly so as to obtain a consistent result. 
A common practice adopted in the literature is to directly deal with Eq.(l), which is a very special 
case of Eq.(3) by simply taking = /ij,"^^* = Q and Ci{jlr) = Ci{Q). Here Q is usually taken as the 
typical spacelike momentum transfer of the process or a value which minimizes the contributions of the 
loop diagrams. As compensation, one varies the value of Q over a certain range, such as [(5/2,2(5], to 
ascertain the renormalization scale uncertainty. This is the simplest scale setting method. It is often 
argued that by setting and varying the renormalization scale in this way, one can estimate contributions 
from higher-order terms; i.e. a change in the renormalization scale will affect how much of a result 
comes from Feynman diagrams without loops, and how much comes from the leftover finite parts of 
loop diagrams. Because of its perturbative nature, it is a common belief that those scheme and scale 
uncertainties will be reduced after finishing a higher-and-higher order calculation. Especially, because 
of the improvement of loop calculation technologies developed in recently years, many high-energy 
processes involving heavy particles have been calculated up to next-to-next-to-leading-order (NNLO) 
or even higher, which greatly improves our theoretical estimations in comparison with the experimental 
data. However, this ad hoc assignment of scale and its range leads to an important systematic error in 
the present theoretical and experimental analysis. Besides the complexity of higher-and-higher order 
calculations, there are many weak points of this conventional scale-setting method: 

1. Following the discussion below Eq.(l), the fixed-order estimation is renormalization scheme de- 
pendent: different choice of renormalization schemes will lead to different theoretical results. In 
addition to the ad hoc dependence on the choice of Q, 

• It is clearly artificial to guess a renormalization scale and to study its uncertainty by simply 
varying /i^ G [(5/2,2(5]. Why is the scale uncertainty estimated only by varying a factor of 
1/2 or 2, and not, say, 10 times of Q ? For example, Ref. [30] argues that after including 
the first and second order corrections to several deep inelastic sum rules which are due to 
heavy flavor contributions, the renormalization scale fir should be taken as /x^ ~ 6.5 m, if 
taking the typical scale Q to be the corresponding heavy quark mass m. The variation of Q 
allows one to estimate some of the contributions from higher-order terms, however, this only 
exposes the {/3j}-dependent non-conformal terms, not the entire perturbative series. 

• Sometimes, there are several choices for the typical momentum transfer of the process, all 
of which, according to the arguments of the conventional scale setting, can be taken as the 
renormalization scale, such as the heavy-quark mass m, the collision energy of the subprocess 
•\/s, etc. Which one provides the correct theoretical estimate? 

2. There are uncancelled large logarithms as well as the "renormalon" terms in higher orders which 
diverge as (n!/3"a"). The renormalon divergence was discovered in 1970s [31, 32, 33]. It has been 
found that those renormalon terms can give sizable contributions to the theoretical estimates, 
such as e+e" annihilation, r decays, deep inelastic scattering, hard processes involving heavy 
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quark, etc.; a detailed discussion on the renormalon problem can be found in the review [34]. As 
a recent example, for the case of H^-boson plus three-jet production at the hadronic colliders, due 
to the renormalization terms and the uncanceled large logs, Ref. [35] shows that a poor choice of 
the scale using the conventional scale setting method can manifest itself as a strong dependence 
on the ratio of next-to-leading-order (NLO) cross section to LO cross section (the so-called K 
factor), which can even predict unreasonable negative NLO QCD differential cross-sections in 
certain kinematical regions. 

3. By taking the Abelian limit A^^ at fixed = CpUg with Cp = {N^ — l)/2Nc, we can 
transform the QCD case effectively to the QED case [29, 36]. A self-consistent scale-setting 
method should be adaptable to both QCD and QED. This fact can be treated as a criterion on 
whether a suggested scale scale setting is correct or not. Conventional scale setting gives wrong 
results when applied to QED processes: As shown above, there is no ambiguity in setting the 
renormalization scale in QED. In the GM-L scheme, the renormalization scale is the virtuality of 
the virtual photon, which naturally sums all vacuum polarization contributions into the coupling. 
There is thus no reason to vary the renormalization scale /i^'^"'" by a factor of 1/2 or 2, since it 
is already the optimized scale. 

4. As more and more data appear, especially because of the running of the high collision energy and 
high luminosity Large Hadronic Collider (LHC), we need more accurate theoretical estimates to 
suit the needs of those forthcoming high precision data. It would be helpful to know to what 
fixed order we can achieve the desired precision. However, the perturbative series does not appear 
to converge, when using conventional scale setting. The conventional scale setting appears as 
a lucky guess, and we have no strict criteria to ensure the perturbative convergence, which is 
especially because of the renormalon terms or uncanceled large logarithms. Taking the top quark 
pair production as an example, it is found that the total cross-section at NNLO level for the 
(gg)-channel, gg — )■ t + t, by taking the conventional renormalization scale choice of rrit, is about 
50% of the NLO cross-section [37, 38]. On the other hand, the experimental result on the ti total 
cross section has been measured with a precision Aau/ca ~ ±7% at the Tevatron [39, 40] and 
~ ±10% at the LHC [41, 42]. Thus, to derive a more precise perturbative estimation, one would 
need to do even higher order calculations, at least at the NNNLO level, which however is not 
expected to be available in the near future. 

In summary, the conventional scale-setting method assigns an arbitrary range and an arbitrary 
systematic error to fixed-order pQCD prediction. One may argue that the correct renormalization scale 
for the fixed-order prediction can be decided by comparing with the experimental data. But this surely 
is process dependent and greatly depresses the predictive power of the pQCD theory. 

For a general fixed-order calculation, what is the correct "physical" scale or optimized scale? To our 
understanding, it should provide a prediction independent of the renormalization scheme and the choice 
of initial scale /x™'*. In fact, this is a criteria of the renormalization group. A pictorial representation 
of what is the optimized renormalization scale is shown in Fig.(l), where the electron-muon elastic 
scattering through one-photon exchange is taken as an illustration. In the GM-L scheme, the optimized 
scale /Xgp^ = t which corresponds to the initial scale- invariant value aem{t)- This optimal scale t is 
independent to the choice of initial scale to; i.e., any choice of to will lead to the same scale t (and then 
same ttem) as shown by Eq.(2). Moreover, by using the proper scale setting method, such as the Brodsky- 
Lepage-Mackenzie (BLM) [20] method and the Principle of Maximum Conformality (PMC) [37, 43, 44, 
45, 46], the prediction is also scheme independent and the argument of the coupling in different schemes 
have the correct displacement. For example, by using the BLM/PMC procedure, one can obtain the 
well-known one-loop displacement between the argument of the coupling in the MS scheme relative to 
the GM-L scheme [47], af^-^{t) = a^{e-^'H). 
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Figure 1: Pictorial representation of the optimized renormalization scale Ho^t- Taking 
electron-muon elastic scattering through one-photon exchange as an example: In the GM-L 
scheme, the optimized scale is /x^p^ = t which corresponds to the scale-invariant value aem{t)- 
As a comparison, the values of a at fixed-orders for different choice of to = under the 
conventional scale setting {i = 1, 2, 3, ■ ■ ■) are shown by thin-and-solid curves. 



As a comparison, the values of ttem at fixed order for different choice of to = = 2, 3, ■ ■ ■) are 
shown by thin-and-solid curves in Fig.(l). The value of strongly depends on the value of (initial) 
scale to under the conventional scale-setting method. Thus, even if a particular choice of to may lead 
to a value of close to aemit) using conventional scale setting, this would not be the correct answer. 
As one includes higher-and-higher orders, the guessed scale will lead to a better estimate; when doing 
the perturbative calculation up to infinite order, any choice of to will lead to the correct value aem(t) as 
required by the RG invariance. However, if one chooses to = t, the complete all-orders result is obtained 
from the onset. 

Does there exist such an optimized renormalization scale for a general high-energy process in non- 
Abelian QCD? If it does exist, how can one set it in a systematic and process-independent way? 

The attempt to solve the renormalization scale and renormalization scheme ambiguities has a long 
history. Many authors have presented their views on how to find such an optimized scale, which have 
served to help clarify the issues involved. In this report, we summarize the principal ideas and results for 
each topic which can be served as a guide to the original literature. We first collect the improvements 
from a general point of view, not following their development in time but follow the sequence of how the 
renormalization scheme and scale questions are understood. Then, we present a detailed introduction on 
the BLM [20] and its underlying principle, the PMC [37, 43, 44, 45, 46, 48]. We shall show that the PMC 
provides the solution for solving the renormalization scale and renormalization scheme ambiguities. 

In Sec. 2 we begin with the RG equation which governs the running (scale) behavior of the QCD 
coupling as(/ir). For convenience, we extend the RG equation also to know the evolution of the renor- 
malization scheme parameters; i.e. the extended RG equations, which was first suggested by Steven- 
son [16, 17, 18, 19], and later improved by Brodsky and Lu [49]. The extended RG equations provide 
a convenient way for estimating both the scheme- and scale- dependence of the QCD predictions for 
a physical process. Any physical observable is independent of the renormalization scale and renormal- 
ization scheme; this is the main property of RG invariance [3, 4, 5, 6, 7]. We utilize the extended RG 
equations for a general discussion on this point. The solution for a special case in which all scheme 
parameters are set to zero, i.e. the 't Hooft scheme [50], is also discussed. The advantage of the 't 
Hooft scheme is that its coupling is scheme-independent and it gives a precise definition for the QCD 
asymptotic scale under a possible renormalization scheme TZ; i.e., the scale for the 't Hooft scheme 
associated with the 7?.-scheme ^qc~d^ [49]. 

As a natural deduction of RG invariance, in Sec. 3 we discuss the self-consistency requirements, such 
as reflexivity, symmetry and transitivity, which must be satisfied by a scale setting method [46, 51, 52]. 
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The transitivity property is especially important for self-consistent scale setting. The fact that the 
renormalization group is called a "group" is mainly because of such transitivity property [3, 4, 5]. 
These self-consistency theoretical requirements can shed light on the reliability of the scale setting 
method suggested in the literature. 

In Sec. 4 we present a brief summary of some typical scale setting methods, such as the Fastest 
Apparent Convergence (FAC) or more strictly the RG-improved effective coupling/charge method [12, 
13, 14, 15], the Principle of Minimum Sensitivity (PMS) [16, 17, 18, 19], the BLM [20] and the PMC [37, 
43, 44, 45]. The FAC and the PMS are designed to improve the perturbative series either by requiring all 
higher-order terms vanish [12, 13, 14, 15] or by forcing the fixed-order series to satisfy the RG invariance 
at the renormalization point [16, 17, 18, 19]. The BLM and the PMC instead improve the perturbative 
series by absorbing only the n^-terms or the {/3j}-terms of the series into the argument of the coupling. 
Thus, these four scale-setting methods have quite different consequences. It has been found that the 
PMS does not satisfy the RG-properties symmetry, reflexivity, and transitivity, so that the relations 
among different observables depend on the choice of the intermediate renormalization scheme [46, 51, 52]. 
Furthermore, the predicted PMS scale for the jet production from e^e~-annihilation does not yield the 
correct physical behavior; it anomalously rises without bound for small jet energy [53, 54]. At present, 
the BLM is widely adopted in the literature and we will present its features observed and developed 
in recent years. The PMC provides the underlying principle for BLM, since it provides a rule to set 
the BLM scales up to all orders, and they are equivalent to each other through the PMC - BLM 
correspondence principle [44]. Thus, all features observed in the BLM-literature are inherited by PMC. 

The main idea of the PMC is that after proper procedures, all non-conformal {/3i}-terms in the 
perturbative expansion are summed into the running coupling so that the remaining terms in the 
perturbative series are identical to that of a conformal theory; i.e. the corresponding theory with 
{Pi} = {0}. The QCD predictions from PMC are then independent of renormalization scheme, because 
the proper displacement of the scales are included. In fact, this can be shown explicitly by considering 
a generalization of the conventional MS'-scheme for dimensional regularization, the T^^-scheme, where 
a further constant 6 from the 1/e poles is subtracted; i.e. j = ^ + ln(47r) — 7^; — 5. The (5-terms in the 
perturbative series will always accompany {/3j}-terms, and thus the elimination of 5-terms is equivalent 
to the elimination of {/9j}-terms. Therefore the PMC estimate can be achieved directly through a proper 
treatment of 5-terms. This leads to a systematic prescription of setting the scales to all-orders, and 
opens the opportunity to start a program for automatically setting the PMC scales [48, 129]. 

It has been found that PMC satisfies all self-consistency conditions. After PMC scale setting, the 
divergent "renormalon" series does not appear in the conformal series; thus as in QED, the scale can be 
unambiguously set by PMC. The scheme independence can be adopted to derive commensurate scale 
relations among different observables and to find the displacements among the effective BLM/PMC 
scales which are derived under different schemes or conventions. The PMC renormalization scale and 
the resulting finite-order PMC prediction are both to high accuracy independent of the choice of the 
initial renormalization scale /i™'*, consistent with RG invariance. Even the residual scale-dependence 
at fixed order due to unknown higher-order {/3j}-terms is substantially suppressed. Since the PMC 
eliminates a serious systematic scale-error in pQCD predictions, it greatly improves the precision of 
tests of the Standard Model (SM) and the sensitivity to new physics at the colliders. It is the main task 
of this report to present a detailed introduction to PMC by including all its developments and useful 
features, its detailed technologies, and its potential phenomenological applications. 

In Sec. 5 we present some applications of PMC, such as the total cross-section or the forward- 
backward asymmetry for the top-quark pair-hadronic production at the NNLO level. In which, we 
show much more subtle points in applying PMC to high energy processes. 

In Sec. 6 we summarize and present an outlook. 
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2 Renormalization Group Equations 



In addition to the purpose of solving the renormahzation scheme and scale dependence of the pQCD 
process, another important goal of a scale setting method is to improve the convergence of the pQCD 
series. A recent review on the development of the QCD coupling is presented in Ref. [55]. The infrared 
behavior of the coupling, in the space-like and time-like regions are discussed using dispersion theory 
and analytic perturbation theory in Refs. [56, 57, 58]. Here we will concentrate on the behavior of the 
coupling and how to deal with its renormalization scheme dependence, based on the RG equation and 
its extended version. 



2.1 Renormalization Group Equation and Its Extended Version 

Predictions for observables in pQCD are expressed in terms of the renormalized coupling whose values 
at any fixed order depends on which renormalization scheme we choose. Conventionally, the scale 
dependence of the coupling is controlled by the /3^-function, 





/3" = /^.£^F^ =-E/3"F^ , (4) 



where the superscript TZ stands for an arbitrary renormalization scheme, such as MS scheme [59], MS 
scheme [47], MOM scheme [60], etc.. The /S^— functions for any MS'-like scheme are the same [61]. 
The various terms in f3]^, f3^, . . ., correspond to one-loop, two-loop, . . . contributions respectively. 
In general, the {Pf^} are scheme-dependent and depend on the quark mass m'j. According to the 
decoupling theorem [62], a quark with mass m'j 3> /i^ can be ignored, and we can usually neglect m'j- 
terms when m'j <C n'^. Then, for every renormalization scale fir, one can divide the quarks into active 
ones with mj = and inactive ones that can be ignored. Within this framework, it is well-known that 
the first two coefficients (3]^i are universal. Hereafter, we simply write them as f3o and (3i. The {f3^^}i>2 
functions for MS'-scheme up to three and four loops can be found in the literature [63, 64, 65, 66]. For 
convenience, we present the results for any simple Lie gauge group with nj fermions and colors [64]: 

11 4 

= y<^A - -^TpUf, (5) 
A = ^Cl-'^CATFUf-ACpTFrif, 



54 ^ 27 ^ " ^ ' 27 ^ ' 9 ^ ^ ^ 9 

^ ^ , /17152 448 \ ^^2^ / 4204 352 \ 424^ , 

^2^^ /7073 656^ \ ^2^2 2/7930 224^ \ 1232^ . 
,ClT.n, ^ f C3) ^CX{^- lic3) + Cmn) (f? - ^C3 

^ A^A V 9 3 ^ A^A V 9 3 

(^abcd^abcd ^ go 704 



where Ca, Na and Cp-, Np are the quadratic Casimir invariants and dimensions of the adjoint and 
fermion representation, respectively, Tp is the trace normalization of the generators of the fermions, ( 
is Riemann zeta function and c^a/f ^'^^ invariant quartic tensors. The expressions for the latter in 
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any simple gauge group can be found in [67]. For the 5't/c'(^)-color-group with fundamental fermions 
the invariants read: 

n n - AT df^'df'"' _ - QN^ + 18 
^^-2' ^^-^i^' ^V:^- 96iV^ ' 

(labcd^abcd iv(iV2 + 6) df^'^df^'^ {N"^ + 36) ^ 



Na 48 ' Na 24 

In particular, for the S't/c'(3)-color-group, we have [64] 



Na = N'- 1, 



/3o = 11-^^/ (9) 
Pi = 102 -yn^ (10) 



^-TFs 2857 5033 325 , 

= nf + n^f (11 

2 18 54 ^ ^ ' 



P^^^ ~ 29243.0 - 6964.30 Uf + 405.089 n] + 1.49931 n) (12) 

We first present the solution of Eq.(4) at the one-loop level, i.e. the solution with only the first term 
keeping in the right-hand-side. The solution is independent of any renormalization scheme TZ and it 
takes the form 




(13) 



where /i™'* stands for an arbitrary initial renormalization scale. This equation governs the one-loop 
behavior of the coupling. It implies that one can first adopt any initial renormalization scale /i™'* to 
measure the coupling; however after we have summed all {/3o}-terms into ^^(/x™^*), its final value will 
not depend on the choice of /x™'*. This is a crucial point, which can be extended to any loops; i.e. if 
summing all types of {/3j^}-terms into the a^-running through the RG equation, the behavior of Q;^(/Xr) 
will be uniquely fixed and is independent of the choice of /i™'*. The final summed result will be more 
accurate. As will be shown later, this fact agrees with the RG invariance and will be a useful guide for 
setting optimal renormalization scales for any fixed-order calculation. On the other hand, if setting 

A|c.^(.rVexp(-^^), (14) 

one can rewrite as{^r) in terms of an overall (universal) scale Aqcd, without any reference to a specific 
initial scale /i™'*. 

Air 

«s(/^r) = 7 (15) 

/3oln(x#^ 

The value of the dimensional scale I^qcd keeps track of the initial parametrization (/i™'*, ^^(/i.J,'^'*)) and 
is universal and scale invariant; its value is not predicted by the theory but must be extracted from the 
measurement of as at a given reference scale. The value of h^qco is commonly believed to be associated 
with the typical hadron size; i.e. to the energy range where confinement effects set in. In effect, Aqcd 
is the scale at which the coupling approximated by Eq.(15) diverges (Landau ghost [68]). 

As suggested by Stevenson [16, 17, 18, 19], it is convenient to use the first two universal coefficients 
f3o and (3i to rescale the coupling and the scale-parameter in Eq.(4); i.e., by rescaling the couphng and 
the scale parameters as [49] 

a = — — a, and = — In/x^, 
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Figure 2: Pictorial representation of the universal coupling a(r, {q}), where r and {q} are 
independent scheme and scale parameters respectively. 



one can express the RG equation (4) into a simpler canonical form 



da^ 



= ^ = -{a'y l + a'' + ^{a'y + 4'{a'y + ■ ■ ■ , (16) 

where cf' = Pf'Po''^ / P\ for z = 2, 3, ■ ■ ■, respectively ^. 

As an extension of the ordinary coupling, one can further define a universal coupling a(r, {q}) to 
include its dependence on both the scale parameter r and the scheme parameters {q}. A pictorial 
representation of the universal coupling a(r, {q}) is shown in Fig. (2). The universal coupling satisfies 
both the scheme and scale evolution equations [16, 17, 18, 19, 49]. 

Following Eq.(16), the scale evolution equation can be rewritten as 

/3(a, {q}) = — = -a^ 1 + a + cga^ + c^a^ + ■ ■ ■ . (17) 

OT L -I 

The scale-equation (17), similar to Eq.(16), can be used to evolve the universal coupling from one scale 
to another. By comparing Eq.(16) with Eq.(17), setting {cj} = {cf^} , there exists a value of r = tt^ 
for which 

a\Tn)=a{TnAcf}). (18) 

This shows that any coupling a^{TTi) can be expressed in terms of a universal coupling a(r, {q}). Notice 
that the evolution equation (17) contains no explicit reference to QCD parameters such as the number 
of colors or the number of flavors. Therefore, aside from its infinite dimensional character, a(r, {q}) is 
just a mathematical function. Truncation of the {ci}-terms (or {/3i}-functions) simply corresponds to 
the evaluation of a(r, {q}) in a subspace where higher-order {cjj-terms are zero. 
The scheme evolution equation is defined as 

/3„(a(r, {c,}), {c,}) = ^a(r, {q}). (19) 



^Another way to rescale the coupling, which is consistent with the large /3o-approximation [69, 70, 71], has also been 
suggested in the literature [72]: i.e. setting = f^o^?' equation (4) changes to 

where df- = /3f 7/3^+^ for i = 1, 2, • • •. 
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The computation of second partial derivative 



(21) 



implies 

g/3„(a(r, {q}),{q}) ^ g/3(a(r, {q}), {q}) 
dr dcn 
which leads to 

/32(a(r, {c.}), {q}) ^ {c.^) \ ^ 

aa(r, {ci}) V {q}), {q}) / 

Finally, we obtain 

/3„(a(r, {ci}), {ci}) = -/3(a(r, {q}), {q}) ^^^^ ^^^y (23) 

where a(0, {0}) = oo and /3(0, {0}) = stand for the boundary conditions. The lower limit of the integral 
has been set to satisfy the boundary condition /3„(a(r, {cj}), {cj}) = 0{a^^^), i.e. a change in c„ can only 
affect terms of order a"""*"^ or higher [16, 17, 18, 19]. The scheme-equation (23) can be used to relate the 
couplings under different schemes by changing {cj}. Equation (23) can be solved perturbatively with the 
help of the scale-equation (17), which can be used to estimate how the uncalculated higher-order terms 
contribute to the final result. An explicit example for this point will be presented in Sec. 4. 3. 6, where 
the value of R{e^e~) at the four-loop level together with its scheme error analysis will be discussed. 



2.2 Solution of the Scale Equation up to Four-Loop Level 

Since any coupling under any renormalization scheme can be related to a universal coupling a(r, {q}), 
the scale-equation (17) can be solved in a conventional way; i.e. the evolution of the universal running 
coupling can be obtained by integrating Eq.(17), which can be rewritten as 

(24) 

where tq = (/3q//3i) ln(/i™^*)2. Here /x™'* stands for an initial renormalization scale. Up to four-loop 
level, it leads to 

L = C + - + \na + (c2-l)a + ^LZl^lH^^ + 0(a3N /25) 
a 2 

where C is an arbitrary integration constant and 

The value of Aqcd can be extracted from a measurement of the QCD coupling at a given reference 
scale or a QCD measure with mass dimensions such as the pion decay constant 
Eq.(25) may be solved iteratively. In fact, up to four loops, it has the form 

„(Mc.»^^H-^ + || + |i.0(±). (2T) 
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Substituting it into Eq.(25), these coefficients kj (j = 1, ■ ■ ■ ,4) can be determined by requiring all the 
terms with ^-^^ (j < 5) vanish. We finally obtain 

a = ^ + ^{C-\nL) + ^[C^ + C + C2-{2C-\nL + l)\nL-l]+^!^c(c^ + ^C + 3c2-2 



I-C3 



3C^ + 5C + 3c9 



3C-lnL + -|lnL 



\nL^ + 0{ — ]. 



(2J 



One will find that the above four-loop solution agrees with Ref. [74] after proper parameter transfor- 
mations. In fact, the integrated RG equation in Ref. [74] takes the form 



L = — 

/3o 



^ + 61 ln(a*) + ih - hlW) + [I - + I) {a*f 



C* 



(29) 



where the following definitions are adopted 



a = — , bi 



TT 



/3o ' 



In 



A2 



(30) 



which are related to our present adopted definitions through the following relations 



/9, = 4^'+^/3;(j = 0,l,---), a=§a* 

Po 



0*2 





(31) 



It is found that by identifying the integration constant C* = ^ (C — \n , the above four-loop solution 
(28) agrees with that of Ref. [74]. 

The universal coupling has a particularly simple form when all the scheme parameters {q} are set to 
zero (the 't Hooft scheme [50]). The 't Hooft scheme is free of higher-order corrections and its running 
coupling a *^ is governed by the simpler RG equation 



'tH 



'tH 



+ In 



'tH 



l + a 



tH 



(32) 



where 



'tH 



PI 
Pi 



In 



^4 



(33) 



and the integration constant C has been absorbed into the asymptotic scale ^qcd ^r convenience. 
It can be solved perturbatively as described above, being a special case of the solution (28). At the 
two-loop level, it however has an analytic solution which can be written as a function of the scale in 
terms of the Lambert function W{z) [75], which is defined through the equation, z = W{z) exp W{z). 

The 't Hooft coupling presents a formal singularity at L *^ = 0; i.e. a'*^ = a(0, {0}) = 00. Inversely, 
it provides a precise definition for the asymptotic scale; i.e., the 't Hooft scale Ag^^,, which is defined 
to be the pole of the coupling in the 't Hooft scheme, a'^^ = a(/3o//3i ln(/i^/(AQ^^)^), {0}). Note that 
since the absorbed integration constant C is arbitrary, the value of ^qcd is not unique, and there are 
infinite number of 't Hooft schemes, differing only by the value of ^qcd- However, under a specific 
renormalization scheme (7^-scheme), its asymptotic scale can be fixed to be the 't Hooft scale associated 
with the 7?.-scheme ^qc~d^ [49], which enters into both 



a(/3oV/3iln(/i.V(A;^cB'')'),{cf} 



and a 



tH 



(/3oV/3iln(/i.V(A'c^cB^)^),{0} 
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Here the word "associated" means we are choosing the particular 't Hooft scheme that shares the same 
't Hooft scale with any given 7^-scheme. In practice, the 't Hooft scale associated with the 7^-scheme 
^QCD^ can be fixed by setting the integration constant to be Ct^. In fact, by taking the same integration 
constant C-ji for both the 't Hooft scheme and the chosen 7?.-scheme, one can obtain a relation between 
^QCD^ and the asymptotic scale Ag^^ for the 7?.-scheme; i.e. 

^Shl^ = exp (l^C^) A^cD- (34) 
As a special case, by choosing Cj^ = Iio.(3q//3i [47, 73], we obtain 

= (I) K'Sd. (35) 

Such a relation is consistent with the observation shown in Refs. [16, 17] and has lately been observed 
in Refs. [44, 49]. The present definition of Aq^j^ is the conventional one, which is associated with the 
choice of C-jjg = Iu/Sq//?! and is originally suggested by Refs. [47, 73]. There are other choices for 
together with the choice of AMf^, [76, 77, 78], which would be helpful in certain cases. 



2.3 Renormalization Group Invariance 

Grunberg has pointed out that any perturbatively calculable physical quantity can be used to define 
an effective coupling, or "effective charge", by incorporating the entire radiative corrections into its 
definition [12, 13, 14, 15]. The effective coupling satisfies the same RG equation as the usual coupling. 
Thus, the running behavior for both the effective coupling and the usual coupling are the same if their 
RG equations are calculated under the same choice of scheme parameters. This idea has been discussed 
in more detail in Refs. [79, 80]. Such an effective coupling can be used as a reference to define the 
renormalization procedure. For example, the effective coupling from the total hadronic cross section 
in e~^e~ annihilation can be defined as [81] 



R„ 



Re 



71 



(36) 



where R^^+^-{Q^) is the Born result and s = is the squared e^e annihilation energy; the effective 
coupling ttgj from the Bjorken sum rule for polarized electro-production can be defined as 



dx 



gA 

9v 



«9i(<5) 



TT 



(37) 



where = —q^ and is the momentum transfer squared. An important suggestion is that all effective 
couplings must satisfy the RG equation [12, 13, 14, 15]. Different schemes or effective couplings will 
differ through the third and higher coefficients of the {/3j^}-functions, which are scheme TZ dependent. 
Thus, any effective coupling can be used as a reference to define the renormalization procedure. 

Physical results should be independent of theoretical conventions. The RG invariance states that a 
physical quantity should be independent of the renormalization scale and renormalization scheme [3, 4, 
5, 6, 7]. Thus it is helpful to use the extended coupling Eq.(18) which contains both the scheme and 
scale parameters for the discussion. 
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2.3.1 Demonstration of Renormalization Group Invariance 



The RG invariance shows that if the effective couphng a{T-ji, {cf-}) corresponds to a physical observable, 
the result from calculating in any scheme should be independent of any other scale ts and any other 
scheme parameters {c^}; i.e. 



dairn, {cf}) 
dcf 



, [scale invariance] 
. [scheme invariance] 



(38) 
(39) 



Demonstration: We provide an intuitive demonstration for the RG invariance from the extended RG 
equations. Given two effective couplings a^Tn, {cf}) and a{Ts, {cf}) defined under two different schemes 
TZ and S, one can expand a^Tn, {cf}) in a power series of a{Ts, {cf}) through a Taylor expansion: 



a{rn, {cf}) 



a{rs + r, {cf + q}) 
da 




(40) 



where f = r-ji — ts, Ci = cf — cf and the subscript S next to the partial derivatives means they are 
evaluated at the point (r^, {cf}). 

The right-hand-side of Eq.(40) can be regrouped according to the different orders of scheme-parameters 
{cj}. After differentiating both side of Eq.(40) over r^, we obtain 



da{rn,{cf}) d^^+^^a{Ts,{cf})f-' 



Or 



(n+l) 



ni 



E 



5("+'^«(r5, {cf }) r 



+ 



[n 



(41) 



where n stands for the highest perturbative order for a fixed-order calculation. It is noted that Eq.(41) 
can be further simplified with the help of RG equations (17,23). If setting n — )■ oo, the right-hand-side of 
Eq.(41) tends to zero, and we obtain the scale-invariance equation (38). This shows that if a{Tn, {cf}) 
corresponds to a physical observable (corresponding to the case of infinite perturbative series, n — > oo), 
it will be independent of any other scale r^. Similarly, doing the first derivative of a(r7^, {cf}) with 
respect to the scheme-parameter c^, one can obtain the scheme- invariance equation (39). 

In other words, if one uses an effective coupling a{Ts, {cf}) under the renormalization scheme S and 
with an initial renormalization scale {ts} to predict the value of another effective coupling a^r-ji, {cf}), 
the RG invariance (38,39) tell us that 

• if we have summed all types of cf -terms (or equivalently the {/9f }-terms) into the effective cou- 
pling, as is the case of an infinite-order calculation, then our final prediction of a(r7^, {cf}) will 
be independent of any choice of initial scale Tg and any renormalization-scheme S. 

• In any case, one needs to set an initial renormalization scale to initiate a calculation, and the 
actual scale may or may not be equal to such initial scale, depending on which scale setting 
method we choose. According to Eq.(41), for a fixed-order estimation (i.e. n ^ oo), there is 
some residual initial-scale dependence. This is reasonable: as shown by Eq.(40), for a fixed-order 
calculation, the }-terms in even higher orders are unknown which however are necessary to 
cancel the scale dependence from the one-lower-order terms. Those unknown-terms provide the 
scale-error source for the fixed-order estimate under the conventional scale setting method. In this 
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method, by varying the scale to be within several times of the typical momentum transfer of the 
process, one can estimate some of the contributions from the higher-order terms, which however 
only exposes the {/3f }-dependent non-conformal terms, not the entire perturbative series. 

If one can find a proper way to sum up all the known- type of {/9f }-terms into the coupling, 
and at the same time effectively suppress the contributions from those unknown-type of {/3f }- 
terms at higher orders, leading to highly convergent perturbative series, such residual initial scale 
dependence can be greatly suppressed. Then, even for a fixed-order calculation, one can eliminate 
the scale error and get the right estimate for a physical observable. 

If setting all the differences of the renormalization scheme parameters to zero, q = = 1, 2, ■ ■ ■), 
Eq.(40) returns to a scale-expansion series for the coupling expanding over itself but specified at 
another scale; i.e. 



a{Tn. |c, I) - a(r5, |c, |) + — r + - — r + - ^-3 r + . (42) 



Using the RG scale-equation (17), the right-hand-side of the above equation can be rewritten as 
perturbative series of a(r5, {cp}), whose coefficient at each order is a {/3j^}-series. This, inversely, 
tells us which {/3/^}-series controls the running coupling at each perturbative order. 

The above discussions are general, which are also suitable for Abelian QED; i.e. by taking the 
Abelian limit A^c — ^ at fixed aem = CpO-s with Cp = {N^ — ^)/'2Nc and Nc the quark's color number, 
we effectively return to the QED case [29, 36]. 

2.3.2 A Combined Evolution of the Coupling in Scheme and Scale 

One can use Eq.(40) together with the scheme and scale evolution equations (17,23) to evolve any 
coupling a{Ts,{cf}), either the usual one or the effective one, "adiabatically" into another coupling 
(^{ttz, {cf'}), not only in scale but also in scheme. Following the idea of Ref. [49], we show how this can 
be achieved. This can be used to relate any two effective couplings. 

First we expand the coupling a^r-ji, {of-}) as a perturbative series of ^(r^, {cf }): 

a{rn, {cf}) = a{Ts, {cf}) + ha'irs, {cf}) + fsa'{Ts, {cf}) + Ua\rs, {cf}) + ■■■. (43) 

This expansion series itself is not accurate if we truncate the series to a fixed perturbative order; i.e., 
if these two schemes TZ and S are quite different and the two scales t-ji and Tg are also quite different, 
then the series might not be convergent. However, it can give us some RG equation improved relations 
for the scheme- and scale- dependent parameters among different schemes. 

From the scheme and scale evolution equations (17,23), up to order C(a^), we have 

= -a\Ts. {cf}) - a\Ts, {cf}) - c^a^Ts, {cf}) + 0{a% 

= 2a%Ts, {cf}) + 5a\rs, {cf}) + 0{a'), 

= -6a\Ts,{cf}) + Oia'), (^^j = a'irs, {cf}) + Oia'), 

= lArs,{cf}) + 0{a'), (-^] =-3a\rsdcf}) + 0{a'). 
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Then, the Taylor expansion of a(r7^, {cf'}) over a{Ts, {cf }) as shown by Eq.(40) can be simphfied as 
airn,{cf}) = a{Ts, {cf}) - fa^Ts, {cf}) + [c2 - f + f^) a\rs, {cf}) + 

'\c, - (cf + 3c2) f + f^] a\Ts, {cf}) + 0{a'). (44) 

After an order-by-order matching of Eq.(43) and Eq.(44), we obtain 

rn = rs-f2, cf = cf - h - + h, cf = Cf - 2hcf + fj + Afl - QM, + 2U. (45) 

Finally, substituting these parameters into the extended RG equations for aij-ji, {cf})-, we can obtain 
the required accurate scheme and scale behaviors of a(r7^, {cf})- 

This finishes the process of deriving the coupling a(r7^, {c^}) at any scale and any scheme from an 
initial or known a(r5, {cf }), which can be summarized in the following [49] : 

• Derive the initial scheme parameters of as {ts, cf, cf , ...) by calculating the coefficients of its 
fundamental f3{as, {cf}) -functions. 

• From Feynman diagram calculation, obtain the expansion series of a{Tn, {cf}) in terms of a{Ts, {cf}), 
e.g. to derive the expansion coefficients /j for Eq.(43). 

• Use the relations as Eq.(45) to identify the parameters tti, cf, cf, ... 

• Use the extended scale evolution equation (17), to run a-n to any required scale. As a byproduct, 
one can adopt the extended scheme evolution equation (23) for an error analysis of a-ji due to 
unknown scheme parameters. 



3 Self-Consistency Conditions for a Scale-Setting Method 

As has been discussed above, the goal of a scale-setting method is to find an optimal renormalization 
scale which can be systematically set in a process independent way. This universality can provide a 
renormalization-scheme independent, initial-scale independent, and also highly convergent perturbative 
series. In the literature, it has been suggested that some self-consistent requirements, such as the 
reflexivity, the symmetry and the transitivity, can shed light on the reliability of the scale setting 
method [46, 51, 52]. These self-consistency requirements have a solid background, which are natural 
requirements of the renormalization group (RG) equation and the RG invariance. 

If one knows how to set the optimal scale, then one can translate the result freely from one scheme to 
another scheme through scale relations [83, 84]. This observation has been emphasized in Ref. [23], where 
the scale transformation among different schemes are called "commensurate scale relations" (GSRs). 
It shows that even though the expansion coefficients could be different under different renormalization 
schemes, after a proper scale setting, one can find a relation between the effective renormalization scales 
which ensures that the total result remains the same under any renormalization scheme. For simplicity, 
following the suggestion of Refs. [46, 51, 52], we also omit the scheme parameters in the coupling in 
the following discussions for the self-consistent requirements of a scale setting method, but will retrieve 
them when necessary. 

In the following four self-consistency requirements are listed, which follow from the RG equation 
and RG invariance: 

1. Existence and uniqueness of the renormalization scale /i^. Any scale setting method must satisfy 
these two requirements. This agrees with our common belief that there does exist an unique 
and optimal renormalization scale for a fixed-order estimation. A pictorial representation of the 
optimized renormalization scale is shown in Fig.(l). For example, the optimal scale for the Abelian 
QED case is set by the GM-L scheme [22]. 
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2. Reflexivity. Given an effective coupling as{fir) specified at a renormalization scale fir, we can 
express it in terms of itself but specified at another renormalization scale /ij., 

a.inr) = + /i(/ir, H , (46) 

where oc ln(/i^//i^^). Up to infinite orders due the scale-invariance (38), we have 

da s{fir)/d In fi^ = 0. This, inversely, means that if as{fir) is known (say, a experimentally 
measured effective coupling), and we try to use the above perturbative equation to "predict" 
as{f^r) from itself, then any deviation of /i^ from /_v would lead to an inaccurate result due to 
the truncation of expansion series. More explicitly, for a fixed-order expansion with the highest 
perturbative-order n, from Eq.(41), we obtain 

OC — 

91n/i;2 n\ 9(ln/i;2 )("+!)■ 

This shows, generally, the right-hand-side of Eq.(46) depends on /^^ at any fixed order. Thus to 
get a correct fixed-order estimate for ^^(/ir), a self-consistency scale setting must take the unique 
value fi'^ = fij. on the right-hand-side of Eq.(46). If a scale setting satisfies this property, we say it 
is reflexive. 

3. Symmetry. Given two different effective couplings asi(/ii) and as2(Ai2) under two different renor- 
malization schemes and at the two renormalization scales fii and fi2 respectively, we can expand 
any one of them in terms of the other: 

asl(/il) = Ois2{f^2) + ruifil, /i2)«s2(At2) H , 

as2{fJ'2) = asiif^i) + r2i{fJ^2, f^i)a1i{f^i) ^ • 

After a general scale setting, we have 

a,i(;Ui) = a,2(/i2) + n2(/ii,/i2)«52(/^2) H > (47) 

as2(yU2) = a,i(/i*) r2i(/i2,/iD"si(/"i) • (48) 

Here as a general choice, we have implicitly set the effective scales at NLO-level to be equal to the 
LO ones; i.e., the effective scales for the highest-order terms are usually taken as the same effective 
scales at the one-lower-order, since they are the scales strictly set by using the known-terms. 

Setting = A2i/ii and /i^ = Ai2/i2, if 

A12A21 = 1 , (49) 

we say that the scale setting is symmetric. 
Explanation: 

If III = A2i/ii and /i^ = \i2l^2^ we obtain 

asi(yWi) = «s2(A2iyUi) + ri2(yUi, A2i/ii)a^2(A2iAti) H (50) 

as2(/i2) = «sl(Al2/i2) + ^21 (/i2, Ai2/i2)a^l(Al2/i2) H • (51) 

As a combination of Eqs. (50,51), we obtain 

asi(Ati) = asi(Ai2A2i/ii) + [ri2(/ii, A2i/ii) + r2i(A2i/ii, Ai2A2i/ii)] asi(Ai2A2i/ii) H . (52) 
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From the reflexivity property, if a scale setting is symmetric, i.e. satisfying Eq.(49), we will obtain 

n2(/il,/i2) +^2l(/i2,AiD = 0, (53) 

and vice versa. This shows that the symmetry property (49) and the relation (53) are mutually 
necessary and sufficient conditions. 

4. Transitivity. Given three effective couplings asi(/ii), as2(/^2); and OssifJ's) under three renormal- 
ization schemes, we can expand any one of them in terms of the other; i.e. 

ttsl(yUl) = as2(/i2) + n2(Aii, /i2)a^2(Ai2) H , 

"s2(yU2) = assifJ^s) + r23{fJ'2, /i3)a^3(/^3) H , 

c^ssifJ's) = o^siifii) + r3i{fi3, fii)ali{fii) -\ . 

After a scale setting, we obtain 



«si(^i) = a,2(/i2) +^i2(/ii,/i2)"^2(/^2) + ■ ■ (54) 

«s2(/i2) = assif^l) + r23if^2, ^l)a1s{^J'l) + ■ ■ ■ , (55) 

«s3(yU3) = a^i(/it) + ri3(/i3,/ii)a^i(/ii) H . (56) 

Setting = A2i/ii, /ig = A32/i2 and /ij = \13H3, if 

A13A32A21 = 1 . (57) 



we say that the scale setting is transitive. 
Explanation: 

If fj,2 = A2i/ii, fJ^l = \z2l^2 and ^\ = Xis^s, we obtain 

asi(/ii) = «s2(A2i/ii) + ri2(/ii, A2i/ii)a^2(-^2i/ii) H , (58) 

as2{l^2) = as3(A32/i2) + ^23(At2, A32/i2)a^3(A32At2) H , (59) 

assifJ^s) = «sl(Al3yU3) +^3l(yU3, Ai3/i3)asl(-^13At3) H • (60) 

As a combination of Eqs. (59,59,60), we obtain 

asl(Atl) = asl(Al3A32A2l/il) + «sl(-^13-^32A2l/il) X 

[r3l(A32A2l/il, Ai3A32A21yUi) + f23(A21yUi, A32A2l/il) + f 12(^1, A21yUi)] H . (61) 

If a scale setting is transitive, i.e. satisfying Eq.(57), we obtain from the reflexivity property, 

ri2(/Ui, /U2) + ^23(/i2> /is) + ^3i(/i3' An) = 0) (62) 

and vice versa. This shows that the transitivity property (57) and the relation (62) are mutually 
necessary and sufficient conditions. The transitivity property shows that under a proper scale 
setting method, we have A21 = A23A31, which means that the scale ratio A21 for any two couplings 
asi and as2 is independent of the choice of a intermediate coupling under any renormalization 
scheme. Thus the relation between any two observables is independent of the choice of renor- 
malization scheme. In fact, the transitivity property provides the theoretical foundation for the 
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Figure 3: Pictorial representation of the self-consistency of the scale setting method through 
the universal coupling a(r, {cj}). The point A with a closed path represents the operation 
of reflexivity. The paths BC and CB represent the operation of symmetry, and the paths 
DF,FE and DE represent the operation of transitivity. 



existence of commensurate scale relations among different physical observables [23]. The transi- 
tivity property is essential for self-consistent scale setting, and is a natural requirement from the 
RG invariance. It has already been pointed out that the transitivity property is the main reason 
why the renormalization group is called a "group" [3, 4, 5]. The transitivity property (57) can be 
extended to an arbitrary number of couplings; i.e. if we have n couplings which are related with 
similar manner as above, then their transitivity relation is 

•^ln'^n(n-l) ' ' ' -^32-^21 = 1- (63) 

One may observe that the Symmetry is a special case of Transitivity, since if setting ass^fJ^s) = 
Q^si(/^i); we have An = 1 and fii(/ii,/ii) = due to the reflexivity, which thus changes the 
transitive relation A13A32A21 = 1 into the symmetric relation A12A21 = 1. 

We present a more intuitive explanation of these requirements based on the universal coupling 
a(r, {cj}) and the extended renormalization group Eqs.(17,23). In the extended RG equations (17,23), 
there is no explicit reference to the QCD parameters, such as the number of colors or the number of 
active-flavors. Therefore, aside from its infinite dimensional character, a(r, {q}) is just a mathematical 
function like, say, Bessel functions or any other special functions [49]. In practice, due to the unknown 
higher order scheme parameters {q}, we need to truncate the beta function /3(a, {q}) and solve the 
universal coupling a(r, {q}) in a finite-dimensional subspace; i.e. we need to evaluate a(r, {q}) in a 
subspace where higher order {cj}-terms are zero. In principle, this function can be computed to arbitrary 
degree of precision, limited only by the truncation of the fundamental /3-function. In this formalism, 
any two effective couplings can be related by some evolution path on the hypersurface defined by 
a(r, {cj}). In Fig. (3) we illustrate the paths which represent the operations of reflexivity, symmetry 
and transitivity. We can pictorially visualize that the evolution paths satisfy all these self-consistency 
properties. A closed path starting and ending at the same point A represents the operation of identity. 
Since the predicted value does not depend on the chosen path, if the effective coupling at A is a a, after 
completing the path we will also end up with an effective coupling a^. Similarly, if we evolve as at B 
to a value ac at C, we are guaranteed that when we evolve ac at C back to the point B, the result will 
be as. Hence, the evolution equations also satisfy symmetry, transitivity follows in a similar manner; 
i.e. going directly from D to E gives the same result as going from D to E through a third point F. 
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In summary, a scale setting method that satisfies uniqueness of the renormahzation scale, reflexivity, 
symmetry, and transitivity effectively establishes equivalent relations among all the effective couplings, 
and thus, among all the physical observables. 

4 Typical Scale-Setting Methods and Their Properties 

According to the RG invariance, physical quantities are renormahzation scheme and scale independent. 
The exact renormahzation scheme independence is respected only approximately for a perturbative 
calculation, which is the well-known renormahzation scheme ambiguity. A resolution of renormahzation 
scheme ambiguity is not simply to find a "good expansion parameter for QCD". In fact, we should 
find a method that can provide the same estimate under any renormahzation scheme for a fixed order 
calculation. There are some suggestions for such purpose, such as FAC [12, 13, 14, 15], PMS [16, 17, 18, 
19], BLM [20] and PMC [37, 43, 44, 45, 48], which are also programmed to solve the renormahzation 
scale ambiguity. A short review of FAC, BLM and PMS can be found in Ref. [85]. Even though 
all of them strive to eliminate the renormahzation scheme ambiguity, they can lead to quite different 
results. For instance FAC and PMS are programmed to directly deal with the nature of the perturbative 
series, whose scale is determined by the total correction; BLM and PMC are programmed to improve 
the behavior of the coupling by absorbing only the part of the correction that is related to coupling 
constant renormahzation (i.e. the rij-terms or {/3j}-terms) into it, which then naturally improves the 
convergence of the perturbative series. 

In this section, we make a detailed discussion on the scale setting methods FAC, PMS, BLM and 
PMC. We present their ideas and basic properties, and show how the self-consistency conditions, such 
as reflexivity, symmetry, and transitivity, are satisfied or broken by these methods. 

4.1 The Fastest Apparent Convergence: FAC Scale-Setting 

It is observed that the standard pQCD prediction for a physical quantity a = /{firf-^QCo) usually gives 
its asymptotic expansion in powers of 1/ In {fir/ Aqcd), which, inversely, means {firf-^QCo) = 
Based on this fact, FAC uses this fact to select the renormahzation scale [12, 13]. The advantage 
of dealing with the inverse function f~^{a) other than the perturbative function f{fir/AQCD) lies in 
that it allows one to get rid of any ambiguity related to the definition of Aqcd and /i^, since any 
redefinition of one of these two scales results only in a trivial overall rescaling of the inverse function 
/~^((t). The inverse function /^^(cr) depends solely on the physical quantity considered, and is therefore 
a renormahzation scheme independent object ^. 

In practice, the FAC scale is determined by requiring all the higher order coefficients in Eq.(l) to be 
zero; i.e. Ci(>i)(/if^^) = 0. It is for this reason, Stevenson called it "Fastest Apparent Convergence" 
(FAC) [17]. It has been argued by Grunberg [14, 15] and Krasnikov [86] that it is really a renormahza- 
tion group improved effective charge or effective coupling scheme, all the known-type of higher order 
corrections can be absorbed into an effective coupling through the RG equation in order to provide 
a reliable estimate, and this method is also applicable when there are large higher order corrections. 
Here, for simplicity, we follow Stevenson's naming for the method. 

4.1.1 Basic Arguments of FAC 

The expansion of /~^(cr) is obtained by introducing an effective coupling asifJr) of the particular 
renormahzation scheme where all higher order corrections to a vanish [12]. If a physical observable in 

^Practically, the inverse function f^^{a) is only an approximation due to a fixed-order calculation, there is residual 
scheme-dependence from the omitted higher order terms. 
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an arbitrary renormalization scheme can be written as 

a = A + B [as{^.r)f 1 + (^l{|^r)o^s{^^r) + 0{al) 
the effective couphng as{iir) is defined by the identity 

a = A + B[as{^^r)]\ 



(64) 



(65) 



where A and B are general perturbative or non-perturbative quantities predicted in principle by QCD, 
d is the Os-order at the Born level and o"i(yUr) is the NLO coefficient. Consequently, as{fir) is the object 
effectively extracted from a LO analysis of the experimental data on a. Next, we require such effective 
coupling also to satisfy the conventional RG equation; i.e. putting p = as/(47r), we have 



2 dp 



/3(p) = -/3oP^-Ap3 + 0(p^). 



Its solution up to two-loop level is [12, 13] 
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For example, we can take Aqcd as Aqcd which is compensated by the MS value of c\^^ to give a 
scheme independent estimate; i.e. 
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(69) 



The result is independent of the choice of MS. Comparison of Eq.(65) with Eq.(67) yields the inverse 
function f~^{a). A simple two- loop approximation is obtained by dropping the integral in Eq.(65), 
giving 

which equals to the two-loop expression in solution (25) under suitable parameter transformations. 
Later on, a more complicated RG equation improved analysis was done by including the three-loop 
/32-term [14, 15]. 

In this way a systematic expansion of {pr/ Aqcd) as a function of p (or equivalently, a) is achieved. 
Since dimensional transmutation is implemented in most direct manner, the only free parameter in 
Eq.(65) is Aqcd (the value of ci is just to compensate the choice of Aqcd in order to provide a scheme- 
independent estimate at the considered perturbative order). Some more points regarding the FAG scale 
setting method are [13]: 

• When the NLO correction in Eq.(68) is large, the use of Eq.(69) amounts to a resummation of 
the most important higher order corrections into Eq.(68): a RG improved perturbation theory is 
achieved. This is the main point of FAG. 

• From the RG equation, assuming the asymptotic expansion of P{p) is well-behaved, an unam- 
biguous criterium for the validity of perturbative theory for each process is given by the condition 
that {Pi/Po)p ^ 1- With the help of Eq.(69), this alternatively means 



Q > Aqcd exp 
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• The FAC method depends sensitively to which quantity it is apphed. For instance, the prediction 
for the ratio i? = (Ti/(T2 of two cross sections 0\ and 02 depends on whether the RG improvement 
is apphed separately to oi and (T2, or directly to R. The first method is more reliable, since a\ 
and o"2 are related directly to Feynman diagrams, whereas i? is a more artificial construct; this 
point of view also has the advantage of exploiting more completely the information contained in 
the expansion of a\ and 02- 

4.1.2 Properties of FAC 

It is straightforward to verify that FAC satisfies all the mentioned self-consistency requirements. 

1. The existence and uniqueness of the renormalization scale fir are guaranteed, since the scale-setting 
conditions for FAC are often linear equations in In/i^, especially for lower order calculations. 

As a simple explanation, if the NLO coefficient Ci{fir) for a physical observable, as defined in 
Eq.(l), has the form 

Ci(/i,) = A + 51n/i^. (71) 
The FAC scale is obtained by requiring, Ciifi^^^) = 0, which leads to 

2. The FAC requires all ln(/i^//i^^)-terms in Eq.(46) to vanish, thus we obtain /i^ = fir- Then, the 
reflexivity is satisfied by FAC. 

3. Symmetry is trivial. After FAC scale setting, two coefficients fi2 and r2i which are defined in 
Eqs. (47,48) satisfy 

n2(/il,/i2) +?'2l(/i2,K) = 0, 

where fil = Ai2/i2 and fi2 = A21/U1. It shows that these two NLO coefficients fi2 and f2i only 
differ by a sign. Thus, requiring one of them to vanish is equivalent to requiring the other one to 
vanish. Furthermore, due to the refiexivity property, one can easily obtain A12A21 = 1. 

4. Transitivity is also satisfied by FAC. In FAC the scales are so chosen such that the NLO term 
vanishes; i.e. after FAC scale setting, Eqs. (54,55) change to 

a,i(/ii) = as2{fJ'l) + 0{al2), (73) 

a.2(yU2) = a.sifil) + 0{al^) (74) 

As a combination, we obtain 

asiifii) = assifil) + 0{a1s) , (75) 

where //* = Xisfis, fi2 = A2i/ii and fil = A32/i2- Notice that this last equation does not contain the 
NLO term. Thus, the relationship between fii and fi^ is still given by the FAC condition (i.e., no 
NLO term), even when we have employed an intermediate scheme. These arguments ensure the 
transitive relation, A31 = A32A21, be satisfied. 
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4.2 The Principle of Minimum Sensitivity: PMS Scale- Setting 



An "unphysical" parameter, such as the renormahzation scale or the renormahzation scheme, means its 
value will not affect the true result of a physical observable. For an all-order calculation, it is true due 
to the RG invariance. However, for a fixed-order calculation, there is a remaining dependence on the 
"unphysical" parameters underlying the conventional scale setting, which depends on the perturbative 
convergence of the process. 

The PMS scale setting is designed to eliminate the renormahzation scheme dependence. Given the 
result in some arbitrary initial renormahzation scheme, the outcome of PMS is suggested to be a unique 
and optimum result, which is scheme independent [16, 17, 18, 19]. It is based on the argument that 
if an estimate has to depend on some "unphysical" parameters, then their values should be chosen in 
order to minimize the sensitivity of the estimate to small variations of these parameters; i.e. the scheme 
and scale must be chosen so as to minimize the sensitivity of the estimation to their small variations. 
It has later been argued, cf. Ref. [87], that the perturbative convergence can also be improved by PMS. 
However, in practice this is not fulfilled. 

More explicitly, the PMS requires the truncated series, i.e. the approximant of a physical observable, 
e.g. pn which is defined in Eq.(l), to satisfy the following RG invariance. 
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where r = \\i{pl/ Is?q(jj~,) and j > 2. Here, we have used the following equation, which is a transformation 
of Eq.(23): 
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The functions /3(as), Os, /?,• and etc. are scheme independent. Here for convenience, we have omitted 
the scheme labels in these equations. 



4.2.1 Basic Arguments of PMS 

Every renormahzation scheme corresponds to a different /3(as)-series, and thus a different (effective) 
coupling. The PMS optimization [16, 17, 18, 19] for the perturbative series can be required in the 
variables that control such a scheme, e.g. the subtraction point pr and the scheme dependent coefficients 
/32, /Ss, . . .. For definiteness, following Ref. [16, 17, 55, 87], we adopt the process R = cr(e+e~ — )■ 
hadrons)/(T(e"'"e~ — )■ p'^p~) for an explanation of PMS. Detailed derivation of the process can be found 
in Refs [55, 87]. For self-consistency, we present their main results here, but will transform their 
notations to agree with our present conventions. 

The quantity Re+e- {s) with an arbitrary choice of pr (for the moment different from the total energy 
s) and in an arbitrary renormahzation scheme takes the form 

where Qf stands for the electric charge of / quark. According to PMS, the quantity Re+e-{,s) should 
be renormahzation scheme and renormahzation scale independent even at the fixed order; i.e., it is 
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stationary for the scale parameter r and the scheme-dependent parameters f3i {i > 2). If we neglect the 
masses of active quarks, we obtain the scale-invariant and scheme-invariant equations: 

h^ + /3(as)T^hRe+e- = 0, (80) 



dr da 



s , 



(3{a,) da'-— ^— i?e+e- = 0, (81) 



\d(3j 'Jo [/3(a')]^^"s 

where j = 2, 3, . . ., r = ln(//^/AQ^^). Similar to Eq.(35), Aqcd is related to the conventional Aqcd 
through the following relation 

Eqs. (80,81) can be used, first to obtain r2, ra, ... in an arbitrary renormalization scheme, when we 
know this quantities in a specific renormalization scheme, and then to make the optimal choice for r, 
/32, /^s, . . .. 

From now on we will use the notation a = as/ir. Replacing (78) in (81), asking that these equations 
for a given fif. are satisfied for an arbitrary value of a, we obtain differential equations for r2, and 
etc.. Restricting to r2, and j = 2, we have 

_ 1 dr2 _ 

dr-i 1 1 ^ dr-i 11 , , 

77^ = 7j/3or2 + -A -77^ = -tt:^ 82 
OT 2 16 0P2 16 po 

Integrating the above equations, we obtain 

where p2 and ps are integration constants independent of r, /32, • • • and are scheme independent. They 
can be calculated, e.g. equating f32, r2, to their expressions (3^^, r^^, r^^ in the MS scheme [88]; 
then, we have 

n - ^MS ^ a ^ 
P2 — — -po m -~ 
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P3 = rr-(rr + V^\ + (84) 
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Note that ps turns out to be independent of s, and r2 has the form 



r2 = -^/3oln + -BoT + r^^^ , (85) 



while rs depends on s and r only through r2. 

Using the 3-loop expression for the {/3j}-functions, we have 
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with 



1 , 4/3o + MA + v^) 



In 



D 4/3o + HPi - VD) 



(87) 



and D = Pi — 4/32/3o • Note that the present complex equation (86) is the strict three-loop solution, as 
a simpler estimation, one can use its perturbative expansion (29) to do the following discussion. 

Let us make the same replacement in Eqs. (80,81) and truncate Eq.(78) at C(a^). By requiring 
Eqs. (80,81) be exactly satisfied, we obtain the following equations 
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Eq.(86) gives a as a function of r, and then, Eqs. (88,89) become equations in r and /32, which determine 
the optimal choice of f{s) and (32{s) for every s. We obtain an optimized running coupling a{s) through 
this way, which together with the optimized values f2(s) and ^3(5), can be used to evaluate the quantity 
Re+e-{s). Some more points regarding the PMS scale setting method are: 

• Because of the scheme-equation (81), we can obtain the coefficients r2, r^, and etc. under any 
renormalization scheme, cf. Eq.(84). This equation is helpful for determining the universal inte- 
gration constants p2, Ps, and etc. by using the results derived from the conventional MS'-scheme. 

• Eqs. (88,89) are general under any renormalization, so the derived formula for the optimized scale 
f(s) is also general, different choice of renormalization scheme will lead to different optimized scale, 
but the final result for R{e^e~) will be the same. So a scheme independent estimate is obtained 
using the PMS. This is the key point of PMS. The optimized scale and hence the optimized 
running coupling can be evaluated numerically [87]. 

• Because of the scheme-independence of the effective PMS scale t(s), one can obtain a relation 
between the effective scales under different renormalization schemes, which could be a commensu- 
rate scale relation as suggested by Ref. [23]. However, according to the above derivation, the PMS 
scale is determined as an overall effective scale for all the considered perturbative contributions, 
so one can not obtain a scale relation as simple as that of BLM or PMC. 

• Following the same way, the PMS method can be extended to higher order approximant. 

• It has been argued that by using PMS, there is a strong correlation between renormalization 
scheme insensitivity and good apparent perturbative convergence [87]. As a naive argument, we 
rewrite Eq.(85) as 
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MS 
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It shows that by using the PMS optimized scale, r2 happens to subtract the /3o-term into the 
coupling, which is similar to the description of BLM and PMC, and then the pQCD convergence 
will be improved at this order in a similar way as that of BLM and PMC. In this sense, PMS 
is consistent with BLM or PMC. However, for even higher order calculations, e.g. for r^, the 
question is much more involved and we have no such simple correspondence. 
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Figure 4: The dependence of the PMS scale parameter T2 as a function the external scale 
parameter ri. 



4.2.2 The Properties of PMS 

Unlike the case of FAC, in general, there are no known theorems that guarantee the existence or the 
uniqueness of the PMS solution. In some processes there may not be a minimum or a maximum. 
Although for practical cases, PMS does provide solutions, and when there are more than one solution 
usually only one of them lies in the physically reasonable region [16, 17, 18, 19], these observations alone 
do not prove that PMS will be trouble-free for new processes. 

To discuss PMS properties in a renormalization scheme-independent way, following the suggestion 
of Refs. [46, 51, 52], we adopt the 't Hooft scheme [50] to define the effective coupling. Under the 't 
Hooft scheme, the RG equation (17) simplifies to 

^ = -0^(1 + a), (90) 

(XT 



r = - + ln(— ^). (91) 



whose solution can be written as 

T = , 

a VI + 0./ 

In the above solution, for convenience, we have redefined r as In ( . /^]j2 ) ; where ^qcd the asymp- 

totic scale under the 't Hooft scheme. 

Given two effective couplings ai and a2 under the 't Hooft scheme, they are related by the pertur- 
bative series 

ai(ri) = a2{T2) + (ra - Ti)al{T2) + ■■■. (92) 
PMS proposes the choice of ji2 (or equivalently, T2) at the stationary point, i.e.: 



doi _ Q _ _d_ 
dr2 dr2 



02(^2) + {t2 - ri)al{T2) 



(93) 



With the help of the above RG equation, we obtain 



1 + «2 = ^ T. (94) 

2(ri - T2) 

In order to express T2 in terms of ri, one must solve the last equation in conjunction with 

^ +logf-^) =r2. (95) 



a2 V i + (22 



In Fig. 4 we present the graphical solution of the PMS scale-parameter T2 as a function of the external 
scale-parameter ri. One may observe two points: 



26 



• T2 > Ti — J. Since T2 7^ ti in any cases, so PMS explicitly violates reflexivity. For a fixed-order 
estimation, when one uses an effective coupling to predict itself, the application of PMS would 
lead to an inaccurate result. 

• In the large momentum region (ri ^ 1), we obtain 02(^-2) — )■ 0, and 

r2^Ti- -. (96) 

Under the same renormalization scheme TZ, we have the same asymptotic parameter ^qcd^ for 
both ai and 02- Here ^qcd^ is the 't Hooft scale associated with the 7?.-scheme, where the word 
"associated" means we are choosing the particular 't Hooft scheme that shares the same 't Hooft 
scale with the 7?.-scheme. Then the relation (96) in terms of /ii and /i2 becomes 

/i2 ^^lexp 1^-^^ . (97) 



More generally, it is found that after PMS scale setting, the scale displacement between any two 
scales /ij and /Xj in the large momentum region is 

A,, = ^ ^ exp f- . (98) 



This would mean that 
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This shows that the PMS does not satisfy the symmetry and transitivity requirements. Let us point out 
that adding the scheme-parameter optimization in PMS does not change any of the above conclusions. It 
only makes the solution much more complicated [87]. The inability of PMS to meet these self-consistent 
requirements resides in that the derivative operations in general do not commute with the operations of 
reflexivity, symmetry and transitivity. This shows the necessity of further careful studies of theoretical 
principles underlying PMS [46]. 

As argued in Sec. 3, any truncated perturbative series will explicitly break RG-invariance (38); i.e. 
Eq.(38) can only be approximately satisfied for any fixed-order estimation. The precision depends on 
to which perturbative order we have calculated, the convergence of the perturbative series, and how 
we set the renormalization scale. As shown by Eq.(76), the PMS requires the truncated series, i.e. the 
approximant of a physical observable, to satisfy the RG-invariance near = /xpMS- This provides the 
underlying reason for why PMS does not satisfy the reflexivity, symmetry and transitivity properties. 
Phenomenological problems of PMS will be discussed in section 4.5. 



4.3 BLM Scale-Setting 

BLM is designed to improve the pQCD estimate by absorbing the {/3j}-terms into the running coupling 
using the n/-terms as a guide. The PMC procedure, which we discuss in detail below, provides a rigorous 
setting for the BLM procedure. Since the invention of BLM by Brodsky-Lepage-Mackenzie [20] in 1983, 
the BLM scale setting method has been widely accepted in the literature for dealing with high energy 
processes, such as the e+e~ hadrons, deep inelastic scattering, heavy meson or baryon productions 
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or decays, the exclusive processes such as pion-photon transition form factors. In 1992, Lepage and 
Mackenzie showed that the apparent failure of QCD lattice perturbative theory to account for Monte 
Carlo measurements of perturbative quantities is a result of choosing the bare lattice coupling constant 
as the expansion parameter [89]. As a solution, they suggested an alternative procedure of BLM for 
determining the effective scale in lattice perturbation theory, which greatly enhances the predictive 
power of lattice perturbation theory. Later on, the reliability/importance of BLM has been emphasized 
in Ref. [90], where an interesting feature for the NLO BFKL Pomeron intercept function u{Q'^,0) 
has been found; i.e. after using BLM scale setting, the intercept function uj{Q^, 0) has a very weak 
dependence on the gluon virtuality in comparison with those derived from the conventional scale 
setting under the MOM scheme and MS scheme. 

In addition, many of the favorable features of BLM have been observed in the literature, which will 
be listed in the following subsections. BLM presents a way to resolve the renormalization scheme-scale 
ambiguity, which results in a new criterion for the convergence of perturbative expansions in QED/QCD 
by unambiguously fixing the perturbative coefficients. In addition to eliminating the renormalizat ion- 
scheme dependence, a better convergence has also been observed because of the absence of renormalons. 
More importantly, the renormalization scale can be determined without computing all higher-order cor- 
rections. Thus, the lower-order or even the LO analysis can be meaningfully compared with experiments. 

BLM scale setting is inspired by QED. As has been discussed in the Introduction, the physical 
quantity within the QED framework can be expanded in perturbative series as 

n 

p„ = Co«f„(/i,) + 5:C,(/i,)«^+^(/i,), ip>0) (101) 

i=l 

where Cq is the tree-level term, Cj stands for the perturbative correction, and p is the power of the 
coupling associated with the tree level. For Abelian theory as QED, since the variation of the effective 
coupling is due to vacuum polarization alone, the BLM method reduces to the standard criterion that 
only vacuum-polarization insertions contribute to the effective coupling. That is, after BLM scale 
setting, we have 

p„ = Co<,(K)+Ci(Ki«?;^'(/^r)+c2(/iri«^;:;'(/^r) + ---, (102) 

where all photon self-energy corrections are absorbed into the couplings by an appropriate (unique) 
choice of effective scales fi*, fi**, ■ ■ ■. In fact, as will be shown by its underlying principle, PMC, the 
BLM scales at different orders are determined by dealing with different {/3j}-functions that will emerge 
in higher order calculation [37, 43, 44, 45]. Since all dependence upon the number of the light-fermion 
flavors {uf) usually enters through the photon self-energy in low orders, both the effective scales //*, 
■ ■ ■, and the low-order coefficients Cj are independent oi Uf. The light-fermion loop corrections serve 
mainly to renormalize the coupling constant aem, as expected. Note that different from the previously 
introduced FAC and PMS methods, each perturbative order will usually have its own scale within BLM; 
there is no reason for running coupling at all orders to have the same scale. In fact, by taking the same 
BLM scales for all orders, serious problems occur, cf. Ref. [91]. 



4.3.1 Basic Arguments of BLM 



The BLM scale can be determined order by order in perturbation theory. We take the LO QCD scale 
setting as an explanation; i.e., to the first order, the physical observable can be re-expanded as [20] 
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where /i^ = /i^"^* stands for an initial renormalization scale, which can be taken as the typical momentum 
transfer of the process. The n f term is due to the quark vacuum polarization and we adopt the familiar 
MS-scheme as an illustration. As will be shown later, by taking any other renormalization scheme, 
one can obtain the same estimate for the physical observable through proper scale displacement [23]. 
It shows that even though the expansion coefficients could be different under different renormalization 
schemes, after BLM scale setting, one can find a relation between the effective renormalization scales 
which ensures that the total result remain the same under different renormalization schemes. 

At the NLO level, all Uf terms should be resummed into the coupling. Using the well-known NLO 
as-running formulae 
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— j and C* = —A + B. (107) 

Both the effective scale fi* and the coefficient are Uf independent. The term 33 A/2 in serves to 
remove those contributions which renormalize the coupling constant into the effective coupling. Some 
more points of BLM scale setting are listed in the following : 

• Using BLM scale setting, one eliminates all n^-terms, so two renormalization schemes that differ 
only by an n/-independent rescaling give identical perturbative expansions in as(/i*). Strictly 
speaking, one absorbs those n/-terms that are related to the renormalization into the coupling 
by using the RG equation (4). Thus, the differences between MS and MS, for example, are 
irrelevant in this approach. Note that for higher-order calculation, there are rij-terms that are 
insensitive to the ultraviolet cutoff and thus have no relation to the /3-function of the coupling, 
such as terms associated with Feynman diagrams with light-by-light quark loops. They should be 
identified and kept separate from the BLM scale setting. 

• The LO BLM scale is determined solely by the parameter A, which exactly comes from quark 
vacuum-polarization insertions. After BLM scale setting, perturbation theory can work well under 
high energy processes, such as e~^e~ — >■ hadrons, deep-inelastic scattering, rjc-decays, heavy {QQ)- 
potential and etc., it has been found that the LO terms in for these processes are by 
themselves quite accurate [20]. 

• Using BLM scale setting, the perturbative expansion will be unchanged in low orders as the 
important momenta vary across a quark threshold, since all vacuum-polarization effects due to a 
new quark are automatically absorbed into the effective coupling. This means, we can use a naive 
LO/NLO as-running with the fixed active fiavor number Uf to do the calculation. In fact, after 
BLM scale setting, the value of can be correctly determined [92]. 

• Reactions with gluon-gluon coupling are more difficult to analyze because of the quark loops 
appear in the higher-order corrections to the gluon-gluon vertex as well as in propagator insertions; 
i.e. it is not easy to separate the divergent part of the vertex from the finite pro cess- dependent 
part in a unique and general fashion. For example, the BLM scale which appears in the three- 
gluon vertex is a function of the virtuality of the three external gluons qf, q^, and q^. It has been 
computed in detail in Ref. [93], where, by taking the subprocess gg ^ g ^ QQ as an example. 
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Figure 5: Typical nj-terms for the electron- muon elastic scattering process at ttg^-order, 
where the solid circles stand for the light-lepton loops. Diagrams (a) and (b) are vacuum 
polarization contributions to the dressed photon propagator which will be absorbed into 
ctemit) as shown by Eq.(2). Diagrams (c) and (d) introduce new type of {/3j}-terms and new 
PMC scales must be introduced. 



the authors show that when the virtualities of the gluons are very different, the energy scale for 
the process should be 

2 2 

f4 oc (108) 

where Ig^inl < kmedl < kmaxi' 9max stauds for the maximal virtuality and etc.. Such scale also 
correctly sets the effective number of quarks (nj) which appear in the /3- function controlling the 
three-gluon vertex renormalization. This example shows that it is critical to set the renormaliza- 
tion scale properly; a prediction based on the guessing scale such as fi^ ~ g^^^ will give misleading 
results. 

• As has been mentioned in the Introduction, there can be residual initial renormalization scale 
dependence due to the unknown higher-order {/3j}-terms. For example, for the simpler QED 
process of the electron-muon elastic scattering through the one-photon exchange only, there is one 
type of {/3j}-terms, which can be conveniently summed up to all orders and its renormalization 
scale can be unambiguously set as the virtuality of the exchanged photon as shown by Eq.(2). 
When two or more skeleton diagrams are involved, more than one types of {/3j}-terms will emerge; 
i.e. Fig.(5c,5d) shows the diagrams with two-photon exchange, and there are two types of {Pi}- 
terms which must be absorbed into two different PMC scales. Because of the unknown higher-order 
corrections for these two types of {/3i}-terms, there is still residual initial scale dependence. 



4.3.2 The Properties of BLM 

It is straightforward to verify that BLM satisfies all the self-consistent requirements outlined in Sec. 3. 

1. The existence and uniqueness of the renormalization scale /x^ are guaranteed, since the scale 
setting conditions for BLM are often linear equations in In/x^. As a simple explanation, if the 
NLO coefficient Ci(/Xr) in Eq.(lOl) has the form 

Ci(/i^) = {a + bnf) + {c + dnf)\n^l, (109) 

with a, b, c and d are constants free of n/, the LO BLM scale can be set as 

ln/i^^^ = -^ + 0(«.), (110) 

where the omitted higher-order a^-terms will be determined by n/-terms at the NLO-level or even 
higher levels. 
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2. Reflexivity is satisfied. The BLM requires all ln(/i^//i'^^)-terms in Eq.(46) to vanish, which are 
proportional to rij-terms, thus we obtain 



3. Symmetry is trivial, because after BLM scale setting, we always have 

That is, those two NLO coefficients only differ by a sign. Thus, requiring one of them to be 
rij-independent is equivalent to requiring the other one also to be n/-independent. This argument 
ensures the symmetric relation, A12A21 = 1, be satisfied after BLM scale setting. 

4. Transitivity is also satisfied by BLM. After BLM scale setting, the two coefficients fi2(/ii,/i2) 
and f23{fi2j f^l) ii^ the following two series 

= as2(/i2) +^i2(/^i,/^2)«s2(/^2) + C'(aL) (111) 

and 



as2{l^*2) = as3{f^;) + f23ift*2,fJ^l)aUl^l) + 0{al,), (112) 

should be independent of n/. After substituting Eq.(112) into Eq.(lll), we obtain 

a.i(/^i) = assif^l) + [^i2(/^i, ^2) + r23{fJ'*2, f^l)] (^Isif^s) + 0{al^) . (113) 

We see that the new NLO coefficient [fi2(/ii, /i2) + f^23{f^2y f^s)] ^ill ^Iso be rij-independent, since 
it is the sum of two n^-independent quantities. These arguments ensure the transitive relation, 
A31 = A32A21, to be satisfied after BLM scale setting. 



4.3.3 Commensurate Scale Relation in QCD 

The BLM prediction is renormalization-scheme independent, which is ensured by the commensurate 
scale relation (CSR) [23], i.e. the specific value of the renormalization scale is rescaled according to 
the choice of the scheme so that the final result is scheme independent. All perturbatively calculable 
observables in QCD, such as the annihilation ratio i?e+e-(Q^); the heavy quark potential, the radiative 
corrections to the Bjorken sum rules and etc., can be related to each other at fixed relative scales. The 
CSR for the observables A and B in terms of their effective couplings {a a and as) takes the following 
form [23] 



. , OIb{Qb) , 2\ 

TT 



(114) 



The ratio of the renormalization scales \a/b = Qa/Qb is so chosen that the coefficient ta/b is inde- 
pendent of the number of flavors Uf. This guarantees that the effective couplings for the observables 
A and B pass through new quark threshold at the same physical scale. The value of \a/b is unique at 
LO, and due to the transitivity of BLM, the relative scales must satisfy the relation 

Aa/b = A^/c Ac/B . (115) 

This ensures that predictions in pQCD are independent of the choice of an intermediate renormalization 
scheme C. In particular, the scale-fixed predictions can be made without reference to theoretically 
constructed renormalization schemes such as MS. 
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Figure 6: Leading order CSRs derived by Ref. [23] for various renormalization schemes that 
are defined through corresponding physical observables. Here av{Q) stands for the effective 
couphng defined from the heavy quark potential V{Q'^) = —4:iTav{Q)/Q'^- 



As a simple explanation of CSR, let us now consider expanding any observable or effective coupling 
aA in terms of ay (corresponding to an arbitrary intermediate renormalization scheme V) up to NLO: 



1 + (CvA + DvArif) 



71 



;il6) 



Note Qa is a formal renormalization scale defined by the physical observable through the effective 
coupling a A- According to BLM scale setting, we must shift the scale Qa in the argument of ay to the 
scale Qv = c^^^^Qa [20], and 



1 + rA/v 



TT 



(117) 



where ta/v = Cva + (33/2)i5vp is the NLO coefficient in the expansion of the observable A in scheme V. 
Thus, the ratio for the two relative scales between the observables A and V, \a/v = Qa/Qv, is fixed by 
the requirement that the coefficient ta/v in the expansion of ay is independent of vacuum polarization 
corrections. Similarly, we can compute another observable or effective coupling as an expansion in 
terms of ay : 

1 + tb/v z ^ ^("^ 



TT 



4) 
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where Qy = Qb/^b/V: and again tb/v niust be independent of vacuum polarization contributions. 

We can now substitute and eliminate av(Qv)! which results in the required LO CSR (114) with 
Qa/Qb = ^A/B = ^A/v/^B/v and ta/b = rA/v — ^b/v- Note also the BLM symmetry property 
^A/B^B/A = 1- Alternatively, we can directly compute the commensurate scale Qa = ^a/bQb by 
requiring ta/b to be n/-independent, which is in agreement with the BLM transitivity. 

The CSR given in Eq.(114) provides a practical way to test QCD: One can compare two observables 
by checking that their effective couplings agree both in normalization and in their scale dependence. 
The ratio of commensurate scales Xa/b is fixed uniquely: it ensures that both observables A and B pass 
through heavy quark thresholds at precisely the same physical point. Calculations are often performed 
most advantageously in MS scheme, but all references to such theoretically-constructed schemes may 
be eliminated when comparisons are made between observables. This also avoids the problem that 
one need not expand observables in terms of couplings which have singular or ill-defined functional 
dependence. 

The intermediate renormalization scheme V, which defines an effective coupling ay, can be taken 
arbitrarily; i.e., in addition to our familiar MS scheme, any perturbatively calculable physical quantity 
can be used to define an effective coupling [12, 13, 14, 15, 79, 80, 81]. In choosing any one of those 
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schemes to predict another observable, the argument of the effective couphng is displaced from its 
(formal) physical value according to CSR, which ensures the final prediction to be independent of 
the renormalization scheme. The relative scale for a number of observables, which are summarized in 
Ref. [23], is indicated in Fig. (6). Here for clarity, we have set av{Q) to be the effective coupling av{Q) 
defined from the heavy quark potential ViQ"^) = —Anav{Q)/Q'^- An essential feature of this scheme 
is that there is no renormalization scale ambiguity, since = —t, the photon/gluon virtuality, is the 
optimized scale in the GM-L scheme. 

Because of CSR, there is no difference of which renormalization scheme one chooses to do the 
calculation. A tricky point is that one sometimes can find a proper scheme which makes the expression 
much more simplified and more convergent. In particular, it has been found that up to light-by-light type 
corrections, all terms involving (3, and tt^ in the relation between the annihilation ratio Re+e- and the 
Bjorken sum rule for polarized electroproduction are automatically absorbed into the renormalization 
scales [23]. Then, the final perturbative series becomes quite simple: 

SgAQ) = s^RiQ*) - s^UQn + s^liQ***), (ii9) 

where a = (3Cir/47r)a, Q*, Q** and Q*** are LO, NLO and NNLO ELM scales accordingly, and the 
two effective couplings a/? and Ogi are defined in Eqs. (36,37). This equation is called the generalized 
Crewther relation. The coefficients in CSR can be identified with those obtained in conformally invariant 
gauge theories as proven by Crewther [94, 95, 96, 97, 98]. 

The CSR between observables can be tested at quite low momentum transfers, even at where pQCD 
expansion would be expected to break down [23]. It is likely that some of the higher twist contributions 
common to the two observables are also correctly represented by CSR. In contrast, expansions of any 
observable in a^jg(Q) must break down at low momentum transfer since aMs(Q) becomes singular at 
Q = Aj^. For example, in the 't Hooft scheme [50] where the higher order = for n = 2,3,..., 
^MsiQ) ^ simple pole aX Q = Aq^^^^^ . The CSR allows one to test QCD without explicit reference 
to renormalization schemes such as MS. It is thus reasonable to expect that the perturbative series will 
be more convergent when one relates finite observables to each other. 

As a summary, the key point of the CSR lies in that the scale displacement between different 
renormalization schemes is unique and does not depend on any intermediate scheme. 

The above discussion on CSR is performed at LO level. In general, such scale relation (115) can be 
extended to any perturbative order. That is, even though the scale values maybe changed from their LO 
values because of the higher-order corrections, the relative relation among the scales must be remained 
unchanged due to the transitivity property of the BLM scale setting. As a demonstration, one needs to 
clarify the following two points : 

• The LO BLM scale itself is a perturbative series with higher perturbative terms coming from a 
higher-order calculation of the physical observable, then we need to show that the scale relation 
(115) is always right for LO BLM scale after including those higher-order terms. 

• We should have similar scale relations for other higher order BLM scales, such as NLO, NNLO 
BLM scales. 

After we demonstrate the first point, then the second point can be recursively demonstrated, since 
according to BLM procedure, if the LO terms have been settled down, we can separate them; while the 
remaining higher-order terms can be regarded as NLO correction to those known terms and then the 
previous LO procedures apply; and so on so forth, we can extend the LO demonstration procedure to 
any perturbative order. In the following, we present a demonstration of how the CSR for the LO BLM 
scales is satisfied even by including perturbative contributions up to NNLO. 

In practice, most physical observable in pQCD are computed in MS scheme, with the running 
coupling fixed at a physical scale of the process. For convenience, we adopt the MS'-scheme as the 
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intermediate scheme. Specifically, up to NNLO, the perturbative series for two effective couplings 
ai{Q)/TT and a2{Q)/T^ which correspond to the two physical observables A and B can be written as 



ai{Q) _ a^{Q) 



+ {Ai+BiUf) 



' \ 2 / \ 3 



TT 



TX 



and 



"2(<5) _ a^{Q) ^ ^ ^^^^^ ( «Ms(<5)~ 



{C2 + D2nf + E2n 



+ 



TT TT \ TX I ■ \ TT I 

On the other hand, the effective coupling ai{Q)/'7i can be directly expressed by a2{Q)/7i; i.e., 
= \- {A12 + Burif) ' ''^ 



(Cl2 + I^12r^/ + ^12rij) 



TT TT ' \ TT I \ TT 

As a combination of Eqs. (120, 121, 122), the coefficients Au, Bi2,Ci2, Du and £'12 read 

A12 = A1-A2, 

B12 = Bi — B2, 

C12 = Ci - C2 - 2(Ai - ^2)^2, 

D12 = D, - D2 - 2{AiB2 + A2B,) + AA2B2, 

E12 = El — E2 — 2{Bi — B2)B2, 



(120) 



(121) 



(122) 



(123) 
(124) 
(125) 
(126) 
(127) 



When taking the NNLO perturbative contributions into account to the physical observable, the LO 
BLM scale will have an NLO term. Following the standard BLM procedure, which will be shown in the 
following subsection, we obtain three LO BLM scales Q^^ , Q^^ and from Eqs. (120,121,122) : 
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(128) 
(129) 
(130) 



B/V 



Q2^/Q and Aa/b 



Q\2 IQ1 obtain the required scale relation 



Aa/B — ^A/v/^B/V — ^A/V^V/B-, 



where we have used the above relations (123-127) for B12 and E12. Here, the second step is due to the 
symmetric relation \b/v^v/b = 1- This finishes our demonstration for the first point. 



4.3.4 An Analytic Extension of MS -Scheme 



As has been discussed in Sec. 2, in conventional M^-scheme, the {/3j}-functions depend on the active 
number of "massless" quarks (tt-/), which is usually a step function of the renormalization scale /i^; i.e., 
the quark masses do not enter into the {/3j}-f unctions since the running coupling is mass independent 
due to the decoupling theorem [62]. 

An important property of BLM scale setting, is that the active number of flavors n / which goes into 
the /3- function can be correctly determined by including the quark mass effect. As has been argued by 
Refs. [56, 57, 58, 92, 99, 100, 101], there are a number of reasons to construct an analytic extension of 
the coupling under M5'-scheme, such as : 
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The comparison of as determined from different experiments and at different momentum scales is 
an essential test of QCD. One source of error is the neglect of quark masses and in the subsequent 
running of from the conventional reference scale, the Z-boson mass. 

Lattice calculation for the heavy quarkonium spectra provides a most precise determination of ag 
at low momentum scales, cf.Refs. [102, 103, 104, 105, 106]. It is important to know how finite 
quark mass effects enter into the running to lower and higher energy scales with as small error 
as possible. 

Finite mass threshold effects in supersymmetric grand unified theories are important when ana- 
lyzing the running and unification of couplings over very large ranges, which has been discussed 
in Refs. [107, 108, 109]. 



An analytic extension of the coupling under the MS'-scheme which incorporates the finite- mass quark 
threshold effects has been suggested in Ref. [92], which is called as MS'-scheme. Such an extension is 
obtained by connecting the coupling under the conventional MS'-scheme to the analytic and physically- 
defined V^-scheme through the GSR derived from BLM; i.e. 

«a7s(v) = ) + ^ + ■••, (131) 

where the commensurate scales Q* and Q** are given by [92] 



Q* = Q exp 
Q** = Qexp 



/105 9^ /103 21 N 

lT28 " iv^ + IT92 + Te^V + ^ I V 



(132) 
(133) 



Here for the l^-scheme, we mean that its effective coupling av{Q) is defined from the heavy quark 
potential V{Q^) = —47rav{Q)/Q'^- There is no renormalization scale ambiguity in ay-scheme, since 

= — t, the photon/gluon virtuality, is the optimized scale in the GM-L scheme which automatically 
sums up all vacuum polarization contribution into the coupling (at high order ay is infrared sensitive, 
so it is difficult to adopt as a standard QGD running coupling). Through a proper way, the l^-scheme 
automatically includes the effects of finite quark masses in the same manner that lepton mass appear 
in Abelian QED [92]. So the extended MS-scheme can also include the quark masses by relating the 
MS'-scheme to l^-scheme through GSR. Note in deriving the LO scale Q*, the ay correction to Q* is 
of less importance for our present analysis, so we do not write it out. 

Taking the logarithmic derivative of GSR given by Eq. (131) with respect to InQ, we can obtain the 
{(]^^'^}-like function, which at LO level gives the following relation between the continuous n^j^ and 
Ufy for various quarks : 




^fMsiQ) = rifAQl =^ I 1 + 7^ I , (134) 

where pi = Q'^/mf [for reference, the quark masses (in GeV) we used are m„ = 0.004, rrid = 0.008, 
rris = 0.200, nic = 1.5, = 4.5 and rrit = 175 ]. Adding all flavors together gives the total n^jjg{Q) 
which is shown in Fig. (7). For reference, the continuous is also compared with the conventional 
procedure of taking n/ to be a step-function at the quark-mass thresholds. The figure shows clearly 
that there are hardly any plateaus at all for the continuous njj^{Q) in between the quark masses. 
Thus, there is really no scale below 1 TeV where n^j^{Q) can be approximated by a constant. We 
also note that if one would use any other scale than the BLM-scale for njrjjg{Q), the result would be to 
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Figure 7: The continuous njjjg{Q) in the analytic extension of the MS'-scheme as a function 
of the physical scale Q [92]. For reference the continuous rif is also compared with the 
conventional procedure of taking rij to be a step-function at the quark-mass thresholds. 



increase the difference between the analytic n j and the standard procedure of using the step-function 
at the quark-mass thresholds. 

Numerically, it is found that taking finite quark mass effects analytically into account in the running, 
rather than using a fixed Uf between thresholds, leads to effects of the order of one percent for the one- 
loop running coupling, with the largest differences occurring near thresholds [92]. These small differences 
are somewhat important for observables that are calculated by neglecting quark masses and could in 
principle turn out to be significant in comparing low and high energy measurements of the strong 
coupling. Moreover, the advantage of the modified scheme a^g^scheme is that it provides an analytic 
interpolation of conventional dimensional regularization expressions by utilizing the mass dependence of 
the physical ay-scheme. In effect, quark thresholds are treated analytically to all orders in mj/Q^; i.e., 
the evolution of our analytically extended coupling in the intermediate regions reflects the actual mass 
dependence of a physical effective charge and the analytic properties of particle production in a physical 
process. Just as in Abelian QED, the mass dependence of the effective potential and the analytically- 
extended MS'-scheme reflects the analyticity of the physical thresholds for particle production in the 
crossed channel. Furthermore, the definiteness of the dependence in the quark masses automatically 

constrains the renormalization scale. Alternatively, one could connect MS'-scheme to another physical 
coupling such as an defined from e~^e~ annihilation. 

By utilizing the BLM scale setting, based on the massless Uf contribution, the analytic extension 
of the MS-scheme correctly absorbs both massless and mass dependent quark contributions from QCD 
diagrams, such as the double bubble diagram, into the running coupling. This gives the opportunity to 
convert a calculation made in the MS'-scheme with massless quarks into an expression which includes 
quark mass corrections from QCD diagrams. In addition one can use this procedure to analytically 
discriminate the dependence of the coupling on time-like and space-like arguments. 

4.3.5 BLM Scale-Setting up to Four-Loop Level 

Based on the main idea of BLM, the method can be extended to higher orders in a systematic way, 
only one should be careful of how to deal with the n/-series at each perturbative order consistently. 
Practically, in doing the BLM extension, the following points must be respected; i.e., 

• All n /--terms, which are associated with the /3-function in the renormalization of the coupling 
constant, must be absorbed into the coupling, while those n/-terms that have no relation to the 
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/3-function should be identified and kept separate. After BLM scale setting, the perturbative 
series for the physical observable becomes a conformal series, all non-conformal terms should be 
absorbed into the effective coupling in a consistent manner. 

• There are always new rij-terms (corresponding to new {/3i}-terms) emerging at each perturbative 
order, so we should introduce new BLM scales at each perturbative order so as to absorb all 
n/-terms into the coupling consistently. There is no strong reason to use a unified effective scale 
for all perturbative orders. In fact, as has already been pointed out in Ref. [91] that if taking only 
one effective BLM scale for the whole perturbative series, one can not fix all unknown parameters 
uniquely even at the NLO level, since one does not have enough constraints to achieve the goal. 

• The BLM scales themselves should be a RG-improved perturbative series [20, 44, 37, 45, 23, 91]. 
The length of the perturbative series for each BLM scale depends on how many new n^-terms (or 
{/3j}-terms) we have from the higher-order calculation and to what perturbative order we have 
performed. 

As a combination of these points, it is interesting to find that the BLM scale setting leads to the 
correct expansion coefficients in the "conformal limit". This inversely shows that the later developed 
PMC scale setting is more essential, since, as will be shown in the following sections, by dealing with the 
perturbative series according to {/3j}-terms, the results are optimized and unique, which also provide 
an unambiguous principle to set the BLM scale up to all perturbative orders. 

In the literature, several ways for extending the BLM method beyond the NLO have been suggested, 
such as the dressed skeleton expansion, the large /3o-expansion, the BLM expansion with an overall 
renormalization scale, the sequential BLM (seBLM), an extension over the sequential BLM (xBLM) 
and etc. [69, 70, 71, 72, 91, 98, 110, 111, 112, 113, 114, 115]. However it can be found that the 
purpose of most of these references is just to eliminate the riy-terms (mostly by introducing an overall 
effective BLM scale for all perturbative orders), but they do not respect all the above listed points 
simultaneously, so even though they do make some improvements in understanding BLM, the criticism 
of BLM made in some of these references are incorrect due to improper understanding/use of BLM. 
Especially, they can not obtain the most important BLM feature that the BLM prediction should be 
independent of the choice of initial renormalization scale. If a method to extend the BLM setting 
up to all orders still depends on the choice of initial renormalization scale, then one will still have a 
(transferred) renormalization scale uncertainty, which is, in principle, not an essential improvement of 
the conventional scale setting method. 

On the other hand, the correct and unambiguous way to the BLM scales up to the two-loop level 
has been suggested by Brodsky and Lu in 1995 [23], and recently such way has been improved up 
to four-loop level [44]. These two references, especially the second reference, can be used as a useful 
guidance for setting BLM scales up to any perturbative order. To be a useful reference and to clarify 
some misunderstandings of BLM in the literature, we present the technical details in the following. 

Generally, separating the n/-terms explicitly at each perturbative order, the pQCD prediction for a 
physical observable p can be rewritten as 



P = 



a:(g) + (Ai + A2n;)<+^(g) + (5i + B^nj + B-in])a'l^\Q) + 
(Ci + C2nf + C^n) + Cin))a':+\Q) + ■■] (135) 



where as{Q) = as{Q)/'i^ and the overall tree-level parameter ro is scale-independent and is free of as{Q)- 
Here, Q stands for the initial renormalization scale, n(> 1) stands for the initial a^-order at the tree 
level. After BLM scale setting, all n/-terms in the perturbative expansion can be summed into the 
running coupling. Here, we shall concentrate on those processes in which all rij-terms are associated 
with the {/3i}-terms. There are n/-terms coming from the Feynman diagrams with the light-by-light 
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Figure 8: The finite electron-loop light-by- light diagrams contributing to the muon's anoma- 
lous magnetic moment [116]. Three more are obtained by reversing the direction of the 
electron loop. 



quark loops which are irrelevant to the ultra-violet cutoff in higher-order processes. However, we should 
remind that there may still be large higher-order corrections not associated with renormalization; an 
important example in QED case is the electron-loop light-by-light contribution to the sixth-order muon 
anomalous moment, its Feynman diagrams are shown by Fig. (8), which is of order (a/vr)^ ln(m^/me) 
and is sizable [116]. 

The BLM scales for the pQCD prediction of p can be determined in a general scheme-independent 
way. The generalization of the BLM procedure to higher order assigns a different renormalization 
scale for each order in the perturbative series, which can be fixed order-by-order. We can shift the 
initial renormalization scale Q into effective ones until we fully absorb those higher-order terms with 
rij-dependence into the running coupling. According to the following steps, one can set the BLM scales 
up to NNLO : 



The first step is to set the effective scale Q* at LO 



<(Q*) + Ai<+i(Q*) + (fii + B2nf)a:^\Q*) + (Ci + Csri; + C:,n))oJl+\Q* 



(136) 



The second step is to set the effective scale Q** at NLO 



-.n+2 



(137) 



The final step is to set the effective scale Q*** at NNLO 



P 



<(Q*) + Ai<+i(Q**) + 5i<+2(Q*-) + Cia:+%Q***) + ■ ■ 



:i38) 



When performing the scale shifts Q Q*, Q* Q** and Q** — )■ Q***, we eliminate the n/-terms 
associated with the {/3i}-terms completely. At the same time, we also have to modify the coefficients, 
since the net changes to the coefficients are proportional to {/3j}-functions. Note that to set the effective 
scale for a"+^, one needs even higher order information and here, as has been discussed previously, a 
sensible choice is Q***, since this is the renormalization scale after shifting the scales up to NNLO. The 
effective renormalization scales up to NNLO can be written as 
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where the effective scales Qq ' are determined so as to ehminate A2nf, B2nf and C2n/-terms com- 
pletely, the parameters x and z are used to eliminate the B^nj and the C^nj terms respectively, and 
the parameter y is used to eliminate the C4nj-term. It is found that 
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+ l)(2n + l)Al - 6n(n + 1)^253 + Qti^C^ 
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3(ri + 2)Bl - 6(n + l)AiC^ 



X = (143) 

(144) 

(145) 
(n + l)AiB2 

If X = 0, or ?/ = 0, or 2; = 0, this shows that there is no new {/3j}-terms that will change the value of 
the effective BLM scale. The exponential form shows that after BLM scale setting, it will not change 
the properties of the initial choice of scale; i.e. its space-like and time-like nature will not be changed. 
The step-by-step coefficients are listed in the following 
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1089(n + 1)A^ + 153nA2 + 66(n + \)A^A2 + (45i - \^mB^)n 
66(n + 1)A^ + 19nA2 + 4(n + 1)^1^2 - 4n(52 + 33^3 

- 40392C4n^ + 143748^2^(3 + + 2n^) + %A2n^{Wx + 35937C4 + 



(146) 
(147) 

(148) 



Co. 



22392^4^^ - 52272^2^(3 + 5n + 2n^) - 24:A2n^{-8C2 + 



64^2^2 _ 

504953^) - 13464A2^n(n2 - 3n - 7) + 72AiA2{l + n)(34A2n - 242^3^ + 
121y42^(3 + 2n)) + 3yl2^r2(2857r2 + 3525i(2 + n) - 95832^3(3 + 2n)) 
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One may observe that even though the perturbative coefficients at the final step Ai, Bi and Ci are 
different, the formulas in deriving them are similar, which shows clearly that the BLM scale setting can 
be done in a recursive way. In deriving the above NNLO formulae, the following equation is implicitly 
adopted, i.e. the value of at any scale Q* can be obtained from its value at the scale Q, 

al{Q) + 0{at). (154) 

For even higher-order corrections, we should use the running behavior derived from Eq.(42) to do 
the calculation, since it shows which {/9j}-terms should be kept in the Os-expansion series. 

All perturbative coefficients Ai, B^, Ci and etc. are renormalization-scheme dependent, so different 
renormalization schemes lead to different BLM scales Q*.**.***; however the final result for p should 
be scheme independent due to GSRs among different observables. Calculating the observable p by its 
corresponding effective coupling and changing to be another effective coupling, starting from Eq.(135) 
and following the same procedures, one can naturally obtain the CSRs up to NNLO. Moreover, by using 
the relations between Q*>**>*** ^nd Q, one can find the needed scale displacement among the effective 
scales which are derived under different schemes or conventions so as to ensure the scheme-independence 
of the observables. For example, from the relation between Q* and Q, one can easily obtain the well- 
known one-loop relation for the couphng [20], af^^{e~^^^Q^) = a^^'~^{Q'^). 

4.3.6 Example of BLM scale setting for Re+e- (Q) O't the Four Loop Level 

Table 1: Coefficients for the perturbative expansion of R^+e-iQ) before and after BLM scale 
setting. 





nf = 3 


% = 4 


nf = 5 


A 


1.6401 


1.5249 


1.4097 


B 


-10.2840 


-11.6857 


-12.8047 


C 


-106.8960 


-92.9124 +2k/15 


-80.0075 +K/33 


A 


0.0849 


0.0849 


0.0849 


B 


-23.2269 


-23.3923 


-23.2645 


C 


82.3440 


82.3440 +2k/15 


82.3440+K/33 
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Measurements of the cross sections for electro-positron annihilation into hadrons provides one of 
the most precise determination of a^. The explicit expression for i?e+e-(Q) up to a^-order under the 
MS'-scheme can be found in Refs. [117, 118, 119]. One finds 



Re+e~{Q)=^Y.^. 



1 + (a*^^(Q)) + A (a^^^(g))' + B [a'''{Q)f + C (a'''{Q) 



(155) 



where 



A = 1.9857 - 0.1152n/, 

B = -6.63694 - 1.20013n/- 0.00518nJ- 1.240r/, 
C = -156.61 + 18.77n/ - 0.7974nJ + 0.0215nJ + Kr/. 

Here rj = Cq) / (s J2q , is the electric charge for the active fiavors. The coefficient k, is yet 
to be determineci, whose contribution will be further suppressed by r], so we set its value to zero in 
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the following numerical calculation ^. The values of A, B and C for n/ = 3, 4 and 5 are presented in 
Table 1. At the present perturbative order, all rij-terms in the above equation should be absorbed into 
as- After applying BLM up to NNLO, we obtain 



1 + (af ^(g*)) + A (af ^(g**))' + B {af^{Q***)f + C (af ^(g*** 



(156) 



where all the coefficients and effective scales can be calculated with the help of the formulae listed in 

the Sec. 4. 3. 5. The coefficients are presented in Table 1, slight differences for B and C with varying nj 
are caused by the charge-dependent parameter rj. 

From the experimental value, r^+^-{?>l.QGeV) = ^Re+e-{^l-QGeV) = 1.0527 ± 0.0050 [122], we 
obtain 

^QCD^ = 412l?°?MeV , A^^ = 359tZMeV. (157) 

With the help of the four-loop coupling (28), we obtain a^'^{Mz) = 0.129l|j:[|?^. It is consistent with 
the values obtained from e+e" collider; i.e., af^{Mz) = 0.13±0.005±0.03 by the CLEO Collaboration 
[123] and aY^{Mz) = 0.1224 ± 0.0039 from the jet shape analysis [124]. One may observe that a 
smaller central value of the world average for a^^{Mz) also results from the measurements of r-decays, 
T-decays, the jet production in the deep-inelastic-scattering processes, and from heavy quarkonia based 
on the unquenched QCD lattice calculations [125]. A larger A;^ leads to a larger af^^{Mz), and vice 
versa. For example, if we set a^^{Mz) to the present world average, we obtain 

A'^\n,=5 = 245l?o MeV , AjYs\nj=5 = 213^1' MeV. 

It is found that after BLM scale setting, the perturbative expansion of Re+e-{Q = 31.6 GeV) becomes 
more convergent. In particular, we find Q* = (0.757 ± 0.008) g which leads to af^{Q*)/af^{Q) = 
1.060 ±0.004. 



4.4 The Principle of Maximum Conformality: PMC Scale-Setting 

Since its invention in 1983, the BLM scale setting has achieved much success in dealing with high 
energy processes. As an extension of BLM scale setting, a program to deal with higher order n^-terms 
associated with renormalization up to NNLO level has been proposed in Ref. [91], which suggests that 
one can expand the effective scale itself as a perturbative series. Later on, an enhanced discussion of 
this suggestion up to NNLO level has been presented in Ref. [23], where the rij-term at the NNLO is 
first identified with /^Q-term and then is absorbed into the running coupling However, BLM in its 
previous form is difficult to be applied to even higher order calculations; i.e. it is not clear how to deal 
with the nj-term, the n^-term, etc. in those higher order corrections in its original version [20]. 

The two most important question for extending the BLM scale setting consistently and unambigu- 
ously up to any perturbative order are: 1) how to deal with the rij-series in the perturbative coefficients 
at each order in an unambiguous way, and what is the underlying principle? 2) how to set the pertur- 
bative series in the BLM scales themselves in a consistent order-by-order manner? It has been found 
that the Principle of Maximum Conformality (PMC) provides the foundations underlying BLM scale 
setting [37, 43, 44, 45, 46, 48, 126, 127]. It inherits all the favorable features of BLM, and it provides the 

Recently, k has been calculated by Refs. [120, 121], which will slightly affect our estimations and its UV-finite 
ny-dependent terms will not affect our BLM treatment. 

"^Strictly speaking, it has been observed that such n^-term together with the n/-tcrm at the same order should be 
rearranged into a proper linear combination of /3i-tcrm and /3g-term. The /3g-term is then absorbed into the LO BLM 
scale and the /3i-term is absorbed into the NLO BLM scale accordingly [44]. 
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Choose renormalization scheme; e.g. of (m™'^) 

i 

Choose fi^™^\ arbitrary initial renormalization scale 

\ 

Identify {Pf} — terms using nj — terms 
through the PMC — BLKl correspondence principle 

\ 

Shift scale of ctg to /i™^ to eliminate {Pf} — terms 

\ 

Conformal Series 




Figure 9: A "flow chart" which illustrates the PMC procedure. 



answer for solving the above two points. Ref. [43] gives one suggestion, where one can use a single global 
PMC scale at LO by proper weighting the separate scales for each skeleton graph, such as t-channel or 
s-channel, to deal with the pQCD cross- sect ion. Later on, the PMC scales were shown to be related 
to the BLM scales in the NNLO analysis through a particular matching of the rij-terms to {/3j}-terms, 
called the PMC-BLM correspondence principle [44]. Recently, taking PMC as a flrst principle a new 
scale setting method was presented, which provides a systematic all-orders method to set the effective 
scales and in a way, which can be readily automatized [48]. 



4.4.1 Basic Arguments of PMC 

The purpose of the running coupling in any gauge theory is to sum all terms involving the ;9^-function. 
Here, we show the scheme dependence in the /3-function explicitly, the superscript TZ stands for an 
arbitrary renormalization scheme. In fact, when the renormalization scales at each perturbative order 
are set properly within PMC, all non-conformal {/Sj^} 7^ terms in a perturbative expansion arising 
from renormalization are summed into the running coupling. The remaining terms in perturbative 
series are then identical to that of a conformal theory; i.e., the theory with = 0. Through this 

treatment, the divergent "renormalon" series of order (a"(/3j^)"'n!) does not appear in the conformal 
series. Thus as in QED, the renormalization scales are determined unambiguously by PMC. 

A "flow chart" which illustrates the PMC procedure is presented in Fig. (9). The PMC provides an 
unambiguous and systematic way to set the optimized renormalization scale up to all orders. We flrst 
arrange all the coefficients, which are usually in n/-power-series at each perturbative order, into 
terms or non-{/3j^}-terms. Then, we absorb all {/3j^}-terms into the running coupling. Note that in 
practice we can directly deal with rij-terms of the coefficients without changing them into {/3j^}-terms, 
and eliminate the n^-terms from the highest power to none also in an order- by-order manner; the results 
are the same due to the PMC-BLM correspondence principle [44]. Different types of {/3j^}-term will be 
absorbed into different PMC scales, and the PMC scales themselves will be a perturbative expansion 
series in a^. After these procedures, all non-conformal {/3j^}-terms in the perturbative expansion are 
summed into the running coupling so that the remaining terms in the perturbative series are identical 
to that of a conformal theory; i.e., the corresponding theory with = {0}. 

As a simple explanation of PMC, for the coefficient Ci(/ir) of the pQCD expansion (101) at the NLO 
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level, we have 

where fir stands for an arbitrary initial renormalization scale, the coefficients Cio (/ir) and Cii{fir) are 
n/-independent, Cio = Cio + f Cn, and Cn = -fCn. The LO PMC scale fj.™^'^^ is then set by the 
condition 

= 0. (159) 

This prescription ensures that, as in QED, vacuum polarization contributions due to the light-fermion 
pairs are absorbed into the running coupling. 

4.4.2 PMC - BLM Correspondence Principle 

A procedure for setting PMC scale at LO has been suggested in Ref. [43], which is adaptable to any 
NLO calculations. Its idea is to distinguish the nonconformal terms from the conformal terms by the 
variation of the cross section with respect to ln(/i™'*)^ (/^r"* stands for some initial renormalization 
scale). Given the analytic form of the hard process amplitude or cross section as a series in ^^(yU™^*), 
one can identify the LO PMC scale through the following way: 

1. There is only one type of /^■'^-function, i.e. /3o, emerges at the NLO level. The variation of the 
cross section with respect to ln(/i™*)^ can be used to distinguish the conformal terms versus the 
nonconformal terms. 

2. The identified nonconformal terms always have the form (3q In [pfj/ (/i™'*)^) where p?^ = pi ■ pj are 
scalar product invariants {i ^ j) which enter the hard subprocess. In practice, these terms can 
be identified as coefficients of rij; i.e., the flavor dependence arising from the light-quark loops 
associated with coupling constant renormalization. 

3. The scale is then shifted from /i™'* to /i^^*" in order to eliminate the non-conformal terms in the 
new perturbative series. Thus, when the scale is correctly set, the coefficients of ^^(yU™*") become 
independent of /3o and ln(/i™*-')^. The series is then identical to that of the conformal theory 
where } = 0. 

At LO, there is only the /3o-term, and the nonconformal terms always have the form of Po ln(/i™*)^, so 
one can determine the nonconformal terms exactly following the above procedures. However, at higher 
orders, the ln(/i™*)^-terms usually appear in a power series as /3o ln(/i™'*)^, /3i ln(/i™*)^, /3o ln^(/U™'*)^, 
etc.. So this method is no longer adaptable to deal with the higher-order corrections, because the 
derivative with respect to a single ln(yU™')^ cannot distinguish all the emerged {/3j^}-terms. 

As noted above, the purpose of the running coupling in any gauge theory is to sum up all terms 
involving the /3 ■'^-function, conversely, one can find all the needed {/3j^}-terms at any relevant order from 
the general expansion (154). This fact should be respected in constructing the perturbative {/3j^}-series 
in both the PMC scales and the perturbative coefficients of a physical observable. Furthermore, using 
this fact and also the known relation between {Pj^} and rif, one can obtain the PMC scales from the 
BLM scale setting method. This is the PMC and BLM correspondence principle [44]. Since 
[i > 2) are scheme-dependent, the PMC and BLM correspondence depends on the renormalization 
scheme beyond the two-loop level ^. 

^It is noted that another suggestion to extend the BLM to all perturbative orders have been given in Refs. [72, 98, 128], 
which is close but different than PMC. At each perturbative order additional {/3]^}-teims are considered. One then 
introduces new free parameters to make the correspondence between n/-terms and {/3j^}-terms. 
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More explicitly, up to NNLO, the physical observable p defined in Eq.(135) can be re-expanded in 
{/3j^}-series as, 



P = ro 



a:{Q) + (A? + A%)a:+\Q) + (5? + B^A + Bl^',)a:+\Q) 

+ C',f3f + CiMi + QVo')ar'(Q)l • (160) 



The results for PMC can be naturally obtained from the BLM scale setting through the following unique 
parameter correspondence; i.e., 

Ai = + (161) 

A2 = —Al (162) 

Bi = Bl + 1025° + 1215° (163) 

B2 = -^(195° + 225°) (164) 

53 = ^5° (165) 
2857 

Ci = C° + ^C° + 1122C° + 1331C° (166) 

C2 = -^(5033C°- 3732^3° -4356C°) (167) 

C3 = ^(325C° + 456^3° + 792C°) (168) 

C4 = -^C° (169) 

which are obtained by comparing the Eqs. (135,160) and the four-loop {/3j^}-terms under the MS scheme 
(TZ = MS), whose expressions have been given in Eqs. (9, 10, 11, 12). 



4.4.3 A Systematic All- Orders Method for PMC Scale-Setting 

Recently another way, i.e. a systematic all-orders method, to set the PMC scales has been suggested 
in Ref. [48]. In comparison to the method of using PMC - BLM correspondence principle, this new 
method has the advantage that the explicit ag-expansion of the effective scales is avoided and the scheme 
independence can be derived without introducing the commensurate scale relation. Moreover, the new 
scale setting processes is easier to automatize. 

The starting point for this method is to introduce a generalization of the conventional schemes used 
in dimensional regularization, where logarithmically divergent integrals are regularized by the following 
transformation of the integration measure: 

J d^p ^ J d'^-^'p , (170) 

where /x is an arbitrary mass scale. Divergences are then separated as 1 /e poles, which can be absorbed 
into redefinitions of the couplings. The choice of subtraction defines the renormalization scheme and can 
be chosen at the theorist's convenience. The arbitrary mass scale becomes the initial renormalization 
scale p}^^^ of the running couplings constants. In the minimal subtraction [MS) scheme one absorbs 
the 1 /e poles appearing in loop integrals which come in powers of 
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where Q is some typical scale and c is the finite part of the integral. The widely used MS-scheme differs 
from the MS'-scheme by an additional subtraction of the term ln(47r) — 7^; together with the 1/e pole. 

One can generalize this by subtracting a constant —5 in addition to the standard subtraction (ln47r — 
7£;) of the MS'-scheme. This amounts to redefining the renormalization scale by an exponential factor; 

i.e. (/i'"*) = (^/i™'*^ exp(ln47r — — 5). In particular, the MS-scheme is recovered for 5 = In 47r — 7^;. 
The (5-subtraction defines an infinite set of renormalization schemes called 5-7?.enormalization (7^^) 
schemes. Moreover, since all TZs schemes are connected by scale-displacements, the /9-function of the 
QCD coupling as in any 7?.5-scheme is the same. In this subsection we use for brevity the notations 
a = as/ 4:71 and (3 as the coupling /3-function in any T^^-scheme. 

Using this generalization, it was shown that a physical observable in any T^^-scheme reads: 



Co + Ciai(/ii) + (C2 + PoCi6i) al{fi2 
+ 



Cs + /3^C^5^ + 2(3oC2S2 + (3'oC^Sf alif^ 



(171) 



where /ij = Qe^'^^, the initial scale is for simplicity set to yU™'* = Q and we defined Ci{Q) = Q. An 
artificial index is introduced on each a and correspondingly on each S to keep track of which coupling 
each 6 term is associated with. The more general expansion with higher tree-level power in a can be 
readily derived [129] and does not change the conclusions and results. 

The expression in Eq.(171) exposes the pattern of {/3j}-terms in the coefficients at each order. Since 
there is nothing special about a particular value of 6, one concludes that some of the coefficients of 
the {/3i}-terms are degenerate; e.g. the coefficient of (3oa{Q)'^ and f3ia{Q)^ can be set equal. Thus, the 
7^5-scheme not only illuminates the {/3i}-pattern, but also exposes a special degeneracy of coefficients 
at different orders. Therefore, for any scheme, the expression for p can be put to the form: 



ro,o + ri,oa(Q) + [r2,o + /3or2,i] a\Q) + rg^o + ftra,! + 2/3or3,i + /3or3,2 a^iQ) 



r^fi + /32r2,i + 2/3ir3,i + -/3i/3o^3,2 + 3/3or4,i + 3/3or4,2 + /3o^4,3 



a\Q) + O(a^) (172) 



where rj^o are the conformal parts of the perturbative coefficients; i.e. = rj^o + C^({A})- In particular, 
it follows that ro,o = Co and ri^o = Ci, while the higher order coefficients Ci>2 are identified with the full 
brackets. 

The artificial indices on Cj and 6i in Eq.(171) reveals how the {/9j}-terms must be absorbed into the 
running coupling. The different 6^8 keep track of the power of the 1/e divergence of the associated 
diagram at each loop order in the following way: the (5^a"-term indicates the term associated to a 
diagram with l/e^~^ divergence for any p. Grouping the different (J^-terms one, recovers in the Nc — ?■ 
Abelian limit the dressed skeleton expansion [29]. Resumming the series according to this expansion 
thus correctly reproduces the QED limit of the observable and matches the conformal series with running 
couplings evaluated at effective scales at each order. 

Using this information from the 5fc-expansion, it can be shown that the order a^{Q) coupling must 
be resummed into the effective coupling a^{Qk), given by: 



1 d l3 

rifia{Qi) = ri^oa{Q)-(3{a)r2,i + -/3{a)—r3^2 + 

2 oa 



r2,oa\Q2) = r2,oa^(g) - 2a(g)/3(a)r3,i + 



,2\n-l^"-+^'^ ' 

(n-1). 



n\ (dln/x^)' 

r4,2 H h '[a)rn+2,n , 



)Rk,n] , (173) 
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which defines the PMC scales Qfc, and where we introduced 



rk,o 




3{k - 1)^^ + {k - l){k - 2)^ 
a oa 



(174) 
(175) 

(176) 



Eq.(173) is systematically derived by replacing the In-' Q\/Q^ by R^^j in the logarithmic expansion of 
of^iQk) up to the highest known ^-coefficient in pQCD. The resummation can be performed iteratively 
using the RG equation for a and leads to the effective scales for an NNNLO prediction: 



, Ql 



R^ 



kA 



/^t\a)Rk,2 + ^t\a)Rk,3 



l + A«(a)i?M + (AL^ 



(i?,,2-i?L) + Af)(a)i?2 



(177) 



The final pQCD prediction for p after setting the PMC scales Qi then reads 

P(Q') = ro,o + ri,oa((5i) + r2,oa'((52) + ^3,00^(^3) + r4,oa^(Q4) + C(a^) 



:i78) 



Here Q4 remains unknown, since it requires the knowledge of rs 1 in the coefficient of . It is noted 
that in contrast to the PMC-BLM correspondence principle, in this method, all effective PMC scales 
are resummed at once, instead of a step-wise process. Moreover, the effective scales naturally become 
functions of the coupling through the /3-function, in principle, to all orders. 

This method systematically sums up all known non-conformal terms, in principle to all-orders, but 
is in practice truncated due to the limited knowledge of the /3-function. It is easy to see that the LO 
values of the effective scales are independent of the initial renormalization scale. This follows since 
taking /i™'' 7^ Q, we must replace -R^^i — > R^^i + ln(5^/(//™'*)^ and thus the LO effective scales read, 
lnQfc,Lo/(/^r"*)^ = Rk,i + lnQV(/^r'''*K wherc p'^'^ cancels and Eq.(177) at LO is recovered. More 
generally the effective scales do not depend on the initial renormalization scale at any order if the 
/3-function is known. In practice, since the /3-function is not known to all orders, there is residual 
renormalization scale dependence, which however is highly suppressed. The effective scales contain all 
the information of the non-conformal parts of the initial pQCD expression for p in Eq.(172), which is 
exactly the purpose of the running coupling. 

In a conformal theory, where = {0}, the 5-dependence vanishes in Eq.(171). Therefore, by 
absorbing all dependence into the running coupling, we obtain a final result independent of the 
initial choice of scale and scheme. It is found that the use of TZs scheme allows us to put this on rigorous 
grounds. From the explicit expression in Eq.(171) it is easy to confirm that 



9p& n, ^^ps 



(179) 



To satisfy scheme-invariance of the physical prediction; i.e. dps/ 86 = 0, we must set the scales such that 
P{a) = 0. This is equivalent to setting the {/3j}-coefficients equal to zero and leads to the conformal 
series. Notice that this holds at any order in perturbation theory and is a theoretical requirement, 
different from the physical fact that the all-orders expression for p, being a physical observable, must 
be renormalization scale and scheme invariant; i.e. dp/d/i™'* = 0. It should be emphasized that this 
is not a fixed point expression for a but is a fully conformal requirement; the /3-function must vanish 
identically. This demonstrates the concept of PMC to any order. 
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4.4.4 Systematic All-Orders PMC Scale Setting for Re+e-{Q) 

We take R^+f,-{Q) as an example to show how to do the PMC scale setting by using the systematic 
all-orders method. The perturbative series matches the generic form of Eq.(172). It can be derived by 
analytically continuing the Adler function, D, into the time-like region [120, 121]: 



1 

2Td 



s+it 



(180) 



with D given by: D^Q"^) = 7(a) — /3(a)^n(Q^, a), where 7 is the anomalous dimension of the vector 
field and 11 is the vacuum polarization function. They can be written as perturbative expansions: 



7(a) = ^na^ , n(Q2, a) = n„(Q2 



:i8r 



n=0 



n=0 



It is then easy to show that to order the perturbative expression for i?e+e- (Q) in terms of 7„ and n„ 
reads (we suppress the Q dependence of Ilj in the following to simplify the notation): 



Rp 



iQ) = lo + liaiQ) + b2-mi]a'iQ) + 



73 + /JiHi + 2/3on2 - (3^ 



+ 



74 + (32li, + 2/3in2 + /^oHs - ^/SoA^ 



3/3, 



2 7^^72 



a\Q) (182) 
a\Q) + 0(a5) . 



This expression has exactly the form of Eq.(172), with the coefficients: 



n,o = 7i > 1), ri^i = (^ > 2), rj,2 
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'li-2 (^ > 3), Ti^s 



i-3 



(z > 4). 



;i83) 



The 7i-coefficients contain iVj-terms in the dimensional regularization schemes, but since they are 
independent of 6 to any order, they are kept fixed in the scale-setting procedure. The expression for the 
coefficient 7^ and Ilj can be found in Refs. [120, 121], while the four-loops /3-function is given in Ref. [64]. 
Now we can set the effective scales Qi, Q2 and up to the NNNLO. By convention, the argument 
of a is space-like; however, the vr^-terms appearing in Re+e~ can be avoided by using a coupling with 
time-like argument, leading to a more convergent series [130]. The last unknown scale can in this case 
be estimated [129]. It turns out that Q4 ^ Q which is the value we have used. The final result with 
five active fiavors reads: 



11' 



-Rp 



iQ) 



1.84 



'«s(<52 



1.00 



11.31 



'as{Q4 



TT 



TT 



(184) 



which is explicitly free of any TZs scheme and scale ambiguities up to strongly suppressed residual ones. 

To find numerical values for the effective scales, we must determine the asymptotic scale, A. From 
the experimental value, re+e-(31.6 GeV) = ^i?e+e-(31.6 GeV) = 1.0527 ± 0.0050 [122], we obtain 



ATcd''' = MeV , A^^Sn 



4191^22 MeV. 



(185) 



The asymptotic scale of TZs can be taken to be the same for any 6, since they share the same /3- 
function. The effective scales are found to be: Qi = 1-3 Q, Q2 = 1.2 Q and Q3 ~ 5.3 Q. These values 
are independent of the initial renormalization scale up to some residual dependence coming from the 
truncated /3-function. 

For completeness, we use our final result to predict the QCD coupling at the scale of the Z-boson 
mass, Mz- as{Mz) = 0.132to;oii. This value is consistent with that of subsection 4.3.6. 



47 



We have checked against the QED case, where Re+e- can be seen as the imaginary part of the QED 
four loop IPI vacuum polarization diagram by the optical theorem, and find in this case nearly complete 
renormalization scale independence of all three scales to NNNLO due to the small value of the coupling. 
Numerically, we obtain for three (lepton) flavors: 

^^QED^'(Q) = 1 + 0.24«e(Qi) - 0.08ae(Q2)' " O.lStteWs)' + 0.05ae(Q4)' , (186) 
where = e^/Air and {^, ^, ^} = {1.1, 0.6, 0.5}. 

4.4.5 Discussion on the Factorization and Renormalization Scale Dependence 

When one applies the PMC scale setting to renormalizable hard subprocesses, the initial and final quark 
and gluon lines are taken to be on-shell so that the amplitude of the hard subprocess is gauge invariant. 
Thus, the application of PMC to hard subprocesses does not involve the factorization scale, and no 
single or double logarithms which involve the factorization scale enter. 
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Figure 10: Diagram for calculating the total cross-section of a general hadroproduction 
process, H1H2 — )■ X, where X stands for any final states. It is obtained from the convolution 
of the partonic subprocess cross-section dij with the parton distribution functions. 



However, for a general hadronic production process such as H1H2 — ?■ X {X stands for any final state), 
which is graphically shown by Fig. (10), we need to be careful of how to set both the renormalization 
scale and the factorization scale consistently. It is important for any correct scale setting method to 
derive the full renormalization and factorization scale dependent terms. The two scales are independent 
quantities. The factorization scale is even necessary even in the theory where the coupling does not 
run. The factorization scale is introduced to identify the domain where the initial and final quarks and 
gluons can be treated as on-shell partons. The total cross-section for a hadronic production process can 
be generally written as 

^ r r dsds 

(TmH2^x = H j j (^1' fif) ^ij{s, M, R), (187) 

where xi = s/S and X2 = s/s. The subprocess cross section dij depends on the renormalization scale 
fir and the factorization scale /i/, with the definitions M = jj^j/Q'^ {Q stands for some typical energy 
scale ^) and R = ^l/fij-. Here S denotes the hadronic center-of-mass energy squared and s = X1X2S is 
the subprocess center-of-mass energy squared. The functions fi/Hi^2i^a, fJ'f) (« = 1,2) are the parton 
distribution functions (PDFs) describing the probability to find a parton of type i with a momentum 
fraction between Xa and x^ + dx^ in the hadron Hi 2- 

^For examples, Q = rriH for Higgs production, Q = rrit for top-quark pair production and etc., so as to eliminate the 
large logs. 
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The factorization scale is the scale entering PDF and the fragmentation functions. It is common 
and simple practice to identify the factorization scale with the renormalization scale, i.e. fif = /i^, and 
then to deal with the process following the same way as that of renormalization scale. The factorization 
scale should be chosen to match the nonpertubative bound-state dynamics with perturbative DGLAP 
evolution [131, 132, 133]. This could be done explicitly using nonperturbative models such as AdS/QCD 
and light-front holography where the light-front wavefunctions of the hadrons are known, a recent 
report of which can be found in Ref. [134]. To fix one's attention on the elimination of renormalization 
dependence, one can fix the factorization scale /i/ to be the value that can eliminate large logs, such as 
setting jif = Q. 

It is important to derive the full renormalization and factorization scale dependence, especially 
those terms from fif ^ fi^ (or R ^ 1), in order to achieve the renormalization scale independence in a 
consistent way. 

As a first step, one can use the fact that the total hadronic cross-section is independent of the 
factorization scale, to derive the first derivative of the subprocess cross-section over the factorization 
scale, /ij^(3"jj(s, M, 1); i.e. to retrieve the factorization scale dependence of the subprocess from 

(Tjj(s, 1, 1) by fixing /i^ = /i/ (or R = 1). For convenience, we rewrite Eq.(187) to a simpler notation 
with the direct-product symbol ®: 



(1^ 



The physical hadronic total cross-section will not depend on fif 

d 



f 



0, 



:i89) 



which leads to the equation 



ijk 



d 



as{fJ'f)Pik ® akj{s, M, 1) + fJ'f^^ijis, M, 1) + as{nf)aik{s, M, 1) O Pkj 

Ofif 



®fj/H2{f^f), 

(190) 



where = Q!s/47r and we have implicitly used the DGLAP evolution equation [131, 132, 133], 



2 d 

/^Z Q..2 



2Ji/Hi 
f 



(191) 



The splitting functions up to three loops can be found in Refs. [131, 132, 133, 135, 136]. Equation (190) 
should hold for an arbitrary fif, therefore, the expression in the square brackets should be identically 
zero for any choice of i and j partons, which inversely yields the following "evolution equation" for the 
subprocess cross-section : 



d 



a,As,MA] 



d 



aij{s,M,l) = -a.(/i/)5I 



Pik ® ^kj{s, M, 1) + aik{s, M, 1) ® P, 



kj 



(192) 



We can solve Eq.(192) through an order-by-order manner using the partonic cross- sect ions (Tij{s, 1, 1) 
as the boundary condition. Initially we can set = f^f = Q {M = R = 1). The splitting function can 
be expanded as 



^-^pS\Q) 



42 



Pl'i(^ 



(193) 
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Figure 11: The scale /i^/yi according to the BLM/PMC, PMS, FAC and the usual 
y^-procedures for the three-jet rate in e~^e~ annihilation, as computed by Kramer and 
Lampe [53, 54]. Notice the strikingly different behavior of the BLM/PMC scale from the 
PMS and FAC scales at low y. The PMS and FAC scales increase at low jet virtuality, which 
is the incorrect physical behavior. 



Then we obtain 
aij{s,M, 1) = aij{s, ^' 1) + I ^j^^ ,,2 



d \ 1 



& 



M=l 



:ain"2 



^■fJ ' M=\ 



(194) 

where Lm = In/ij/Q^. 

As a second step, one can adopt the RG equation to retrieve aij{s, M, R) from aij{s,M, 1). This 
step can be done by directly replacing the coupling at the renormalization scale /i/ to /ir; i-e. using the 
following formula. 



(195) 



where = as/n. The full renormalization and factorization scale dependence for the splitting function 
can be retrieved by using this equation; i.e. [136], 



43 



f^r f^r f'r 



0(1)/ 



■(196) 



One can then apply e.g. PMC to set the renormalization scale for the splitting function. 

The above two procedures can be extended up to any perturbative order. After doing these two steps, 
one can obtain the required aij{s, M, R) with full renormalization and factorization scale dependence. 



4.5 A Comparison of FAC, PMS and BLM/PMC 

As shown above, because the scale setting methods, such as FAC, PMS and BLM/PMC, have quite 
different starting points, they can give strikingly different results in practical apphcations. For example. 
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Kramer and Lampe have analyzed the apphcation of the FAC, PMS, and BLM/PMC methods for the 
prediction of jet production fractions in e~^e~ annihilation in pQCD [53, 54]. Usually, jets are defined by 
clustering particles with invariant mass less than ^/ys, where y stands for the resolution parameter and 
^/s is the total center-of-mass energy. Physically, one expects the renormalization scale /i^ to reflect the 
invariant mass of the jets, that is, fir should be of order ^/ys. For example, in the analogous problem 
in QED, the maximum virtuality of the photon jet which sets the argument of the running coupling 
cannot be larger than y/ys. Thus one expects /i^ to decrease as ?/ — )■ 0. However, as shown by Fig. (11), 
the scales chosen by the FAC and PMS methods do not reproduce this physical behavior: The predicted 
scales ^^'^ and fi^^'" rise without bound at small values for the jet fraction y. This shows that the 
FAC and PMS can not get the right physical behavior in this limit, since they have summed physics 
into the running coupling not associated with renormalization. On the other hand, the BLM/PMC 
scale has the correct physical behavior as ?/ — )■ 0. Since the argument of the running coupling becomes 
small using the BLM/PMC method, standard pQCD theory in will not be convergent in the low y 
domain [137]. In contrast, the scales chosen by PMS and FAC give no sign that the perturbative results 
break down in the soft region. 

The real power of FAC is the concept of the effective charge, which allows one to define a scheme 
defined from a physical observable. As we described, the commensurate scale relations can be used to 
relate one effective charge to another. Furthermore, we list the main differences for PMS and BLM/PMC 
in the following: 

• The PMS chooses the renormalization scale such that the first derivative of the fixed-order cal- 
culation with respect to the scale vanishes. However, this criterion of minimal sensitivity gives 
predictions which are not the same as the conformal prediction. As shown in Ref. [87], by using 
the PMS together with the scheme-equations (23) and the scheme- independent equation (39), the 
renormalization scheme dependence can be reduced to a certain degree through an order-by-order 
manner. But there are still residual scheme dependence due to unknown higher order corrections, 
and in principle, the PMS prediction depends on the choice of renormalization scheme ^, and it 
disagrees with QED scale setting in the Abelian limit. Most important, the PMS does not satisfy 
the RG-properties symmetry, reflexivity, and transitivity, so that relations between observables 
depend on the choice of the intermediate renormalization scheme. Hence, when we successively 
express one effective charge in terms of others, PMS would lead to inconsistent scale choices. 

• At the present the BLM has been widely accepted for dealing with high energy processes. The 
PMC, being the principle underlying BLM, sums all {/3j}-terms in the fixed-order prediction into 
the running coupling, leaving the conformal series. The PMC is equivalent to BLM through 
the PMC-BLM correspondence principle. It satisfies all the RG-properties, such as reflexivity, 
symmetry, and transitivity. The PMC prediction is thus scheme-independent, and it automatically 
assigns the correct displacement of the intrinsic scales between schemes. The variation of the 
prediction away from the PMC scale exposes the non-zero dependent terms. The PMC 
prediction does have small residual dependence on the initial choice of scale due to the truncated 
unknown higher order {/3i}-terms, which will be highly suppressed by proper choice of PMC scales. 

5 Applications of PMC 

In this section, we present some recent examples for the PMC scale setting. Some more subtle points 
in using the PMC scale setting are presented, which are useful references for future applications. 

^It is hard to estimate the contributions from those unknown higher-order terms within the framework of PMS. While 
for BLM/PMC, such dependence can be analyzed by using the extended renormalization group [49]. 
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5.1 Top- Quark Pair Total Cross Section at the NNLO Level 

The total top-pair production cross-section au has been measured at the Tevatron with a precision 
Aatt/au ~ ±7% [39, 40] and the two LHC experiments have reached similar sensitivity [41, 42]. The- 
oretically, the total cross-section for the top-pair production has been calculated up to NLO within 
the MS'-scheme in Refs. [138, 139, 140, 141, 142, 143]. Large logarithmic corrections associated with 
the soft gluon emission have been investigated and resummed up to next-to-next-to-leading-logarithmic 
corrections [140, 143, 144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154]. These results provide the 
foundation for deriving a more precise estimation by using PMC. 

5.1.1 Basic Formulas 

According to Eq.(187), The total top-pair production cross-section can be written as: 

^ If f ds 

(TH^H2^tt+X = Y1^ J J —fi/Hi {Xl,fif) fj/H2 {^2, f^f) ^ij{s, (/ir) , /i/) , (197) 

where Xi = s/ S , X2 = s/ s. The top-quark mass rrit is the mass renormalized in the on-shell (pole-mass) 
scheme. Setting s = AmflS/Amfy^ and s = s(S'/s)^^, we can transform the two-dimensional integration 
over s and s into an integration over two variables 2/1.2 G [0, 1], which can be calculated by using the 
improved VEGAS program [155, 156, 157]. The partonic cross-section atj can be decomposed in terms 
of dimensionless scaling functions fij. Up to NNLO, it takes the form 



(Tij 



where p = imf/s, as{Q) = as{Q)/n and (ij) = {{qq), (gg), (gq), (gq)} stands for the four production 
channels respectively. Q stands for the typical energy scale of the process. When applying PMC to 
the renormalizable hard subprocesses which enter the pQCD leading-twist factorization procedure, the 
initial and final quark and gluon lines are taken to be on-shell so that the hard subprocess amplitude is 
gauge invariant. Thus, the application of PMC to hard subprocesses does not involve the factorization 
scale. It is convenient to fix the factorization scale pf = Q. In principle this uncertainty can be removed 
if one knows the bound state wave function. As for the initial renormalization scale /i™'*, we also set 
its value to be Q. 



The scale functions fij^'^{p,Q) can be directly read from the HATHOR program [152]. According 
to the PMC scale setting, we need to find the explicit terms that are nj- or nj- dependent, which 
should be absorbed into the a^- running. The QCD Coulomb-type correction may also provide sizable 
contributions in the threshold region [158, 159], so terms that are proportional to tt/v or {tt/vY {v = 
y/l — p, the heavy quark velocity) should be treated separately [160]; i.e., we need to introduce new 
PMC scales for the Coulomb type terms. Then, the NLO and NNLO scaling functions can be rearranged 
as 

'it 

V 



fl,{p,Q) = [Auj + B,,,nj] + Duj(-) (199) 
4.(p,Q) = \A2ij + B2^Jnf + C2ijn} + [D2^J + E2ijn f] ( - ) + F2^j ( - ) . (200) 



Substituting them into Eq.(198), the partonic cross-section aij changes to 



aij 



AoijaliQ) + l^[Auj + BujUf] + Duj } a^(Q) 



A2ij + B2ijnf + C2ijnj\ + [D2ij + ^2ijri/] f - ) + F2ij ( - ) i a^(Q), (201) 
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where Aoij = f^j{p,Q). 

The PMC scales can be determined in a general scheme-independent way as described in Sec. 4. 3. 5. 
We shift the renormalization scale Q to absorb the {/3j}-dependent terms, using the n/-dependence as 
a guide: 



m^ (Jij 



Aiij 



TX 



D 



lij 



2k 



1 — exp(— 2k) 



«s(Q;),(202) 



where k = ^a,(g^) + §^ (A asiQ*2). Ql'** are LO and NLO PMC scales for the usual part and Ql 
is the LO PMC scale for the Coulomb part. For the usual part, the two PMC scales Ql'** satisfy 
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(203) 
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where the coefficients are [37] 

2Aiij + 335 



Aiij 

A2ij 
B2ij 



lij 



A 



2A 



2ij 



2ij 



335 



2ij 



SA^ij 
1 



4A 



Oij 



[mBuj{2Auj + 335i,,) + Ao^j{%A2^j + 3065i,j - 2178^2,,)] 
[AAQ,,{B2ij + 33C2y) - Bu.ilQAoij + QAuj + ^Buj)] . 



(204) 
(205) 

(206) 



While for the Coulomb part, we have adopted the Sommerfeld-Gamow-Sakharov rescattering for- 
mula [161, 162, 163] to sum up the higher-order vr/w terms. The overall factor (vr/f) before Duj 
shall be canceled by a f ^-factor from the phase space. Its LO PMC scale Q2 satisfies 



In 



Qf_ 
Q' 



2E- 



2ij 
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lij 



(207) 



and the coefficient 

b^ij = {2D2i, + 3352ij)/2. (208) 

Since the channels [ij) = {{qq), [gg), (gq), {gq}} are distinct and non- interfering, their PMC scales 
should be set separately. 

• For (gg) channel, all the coefficients A^gg, Aiqq and etc. are non-zero. Following PMC scale setting, 
the rij-terms that are associated with the {Pf^^}-ieims are absorbed into the Os-coupling: 

— When performing the scale shift Q ^ Ql, the first type of {/3j^}-terms Biqq and C2qq are 
eliminated. Part of i?2qq which contains the same type {/3j^}-term is also absorbed into 
running. The remaining part of -8299 is compensated by Aiqq and Bigg to ensure that the first 
type of {/3j^}-terms are absorbed into -coupling exactly, which results in a new variable 
B2qq. Because [Biqq/ A^qq) shows a monotone increase with the increment of collision energy 
i/s, Ql shall show the same trend versus y/s. 

— When performing the scale shift Ql Ql* , the second type of {/3i}-terms, i.e. B2qq are 
eliminated. As shown by Fig. (12), the value of B2qq is always negative and Aiqq has a 
minimum value in lower a/s, then one can find that the NLO PMC scale Ql* shall be 
suppressed to a certain degree in comparison to Q. 

— The Coulomb type correction provides a distinct contribution to the total cross-section in the 
threshold region, which should be treated separately from the usual part. Similarly, when 
performing the scale shift Q ^ Q2 for the Coulomb type contribution, i?2gg is eliminated. 
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Figure 12: PMC coefficients of tlie [qq) cliannel versus tlie subprocess collision energy ^/s. 
mt = 172.9 GeV. 



For (gg) channel, we have C2gg = 0, while other coefficients AQgg, Aigg and etc. are non-zero. It 
can be treated in a similar way as the {qq) channel. It is found that in distinction to the q + q 
case, both 2B2gg/ Aigg and SBigg/AQgg are close to zero, and thus its LO and NLO PMC scales 
{Ql and Ql*) become close to Q with appropriate modifications due to ^/s. 

For (gq) or (gq) channel, we have = 0, Digg = 0, E2gq = F2gq = 0, which shows that the 
Coulomb type corrections start only at the NNLO order. We only need to set one LO PMC scale 
for these two channels. That is. 



m^aij 



[Auj] aliQ) + [[A2^j + B2.jnf] + [D2.,] } at{Q) 



AujaiiQl) + A2i, atiQl) + D2., - a',{Q) 



TT 



V 



where {ij) = (gq) or (gq), igij = ^2ii + 3352^^/2 and \nQf/Q^ = 232,^/ A 



lij- 



5.1.2 Numerical Analysis for the Total Cross Section 

To do the numerical calculation, we adopt nit = 172.9 ± 1.1 GeV [28] and the CTEQ PDFs of version 
2010, i.e. CTIO [164]. The combined PDF and uncertainty are set by using different PDF sets 
determined by varying as{mz) G [0.113, 0.230]. At present, to keep our attention on the renormalization 
scale, we set as usual, fif = rrit. As initial choice, we set /x™'* = nit. After PMC scale setting, the 
PMC scales are usually different from mt, so one must use the formulas listed in Sec. 4. 4. 5 to get 
the full renormalization and factorization scale dependence before applying PMC. This point is very 
important for the later initial-scale-independent analysis. In the literature, the full renormalization 
and factorization scale dependence for the top-quark pair production up to NNLO can be found in 
Refs. [147, 152]. 

We first present the total cross-sections for the top-quark pair production using the PMC scale 
setting by fixing all the input parameters to their central values. The results are presented in Tables 2 
and 3, where for comparison, the total cross-sections for the conventional scale setting method {fir = ^t) 
are also presented. 

• Tables 2 and 3 show that the pQCD convergence is improved after PMC scale setting. This is 
due to the fact that we have resummed the universal and gauge invariant higher-order corrections 
associated with the {[3^^^}-ieims, into the LO and NLO -terms by suitable choice of PMC scales. 
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Table 2: Total cross-sections for the top-pair production before and after PMC scale setting 



for = 1.96 TeV. rrit = 172.9 GeV and the central CTIO as the PDF [164], 





before PMC scale setting 


after PMC scale setting 




T n 


NT n 


NNT n 

IN IN ijV^ 


Loiai 


T n 


NT n 


NNT n 


LuLd.1 


q -\- q [po) 


4.989 


0.975 


0.489 


6.453 


4.841 


1.756 


-0.063 


6.489 


9 + 9 (pb) 


0.522 


0.425 


0.155 


1.102 


0.520 


0.506 


0.148 


1.200 


g + q (pb) 


0.000 


-0.0366 


0.0050 


-0.0316 


0.000 


-0.0367 


0.0050 


-0.0315 


g + q (pb) 


0.000 


-0.0367 


0.0050 


-0.0315 


0.000 


-0.0366 


0.0050 


-0.0316 


sum (pb) 


5.511 


1.326 


0.654 


7.492 


5.361 


2.188 


0.095 


7.626 



Table 3: Total cross-sections for the top-pair production before and after PMC scale setting 



for ^/i = 14 TeV. mt = 172.9 GeV and the central CTIO as the PDF [164]. 





before PMC scale setting 


after PMC scale setting 




LO 


NLO 


NNLO 


total 


LO 


NLO 


NNLO 


total 


q + q (pb) 


73.45 


9.00 


5.16 


87.58 


69.33 


21.08 


-2.98 


86.96 


9 + 9 (pb) 


487.52 


262.96 


49.87 


800.68 


485.51 


303.69 


37.35 


835.41 


g + q (pb) 


0.00 


9.30 


7.69 


17.00 


0.00 


9.37 


7.69 


16.97 


g + q (pb) 


0.00 


0.02 


1.92 


1.95 


0.00 


-0.03 


1.92 


17.06 


sum (pb) 


560.96 


281.29 


64.63 


907.20 


554.84 


334.11 


43.98 


941.26 



It is also the reason why after PMC scale setting, the total cross-section au is increased by ~ 2% 
at the Tevatron and ~ 4% at the LHC. This small change in the total cross-section after PMC 
scale setting means that the naive choice of fir = rrit is a. viable approximation for estimating the 
total cross-section. 

• Since different channels have quite different behaviors, it is necessary to use different PMC scales 
for each channel. The PMC scales are functions of ^/s, whose behaviors up to ^/s = 7 TeV are 
presented in Fig. (13). Because of the behaviors of the PMC coefficients, the LO PMC scale for the 
(gg)-channel increases with ^/s and is much larger than for large a/s. As a result, its LO cross- 
sections at the Tevatron and LHC are decreased by 3% — 5% relative to the standard guess under 
the conventional scale setting. Because \Bigg/AQgg\ <^ 1, the LO PMC scale for the (5f5f)-channel 




Figure 13: PMC scales versus the sub-process collision energy y/s for the top-quark pair 
production up to y/s = 7 TeV, where we have set the initial renormalization scale /x™'* = Q. 
Here Q = mt = 172.9 GeV. 
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Figure 14: The ratio -R™*" = ^pmO| ^ ^'^ versus the sub-process colhsion energy -y/i up to 

14 TeV, where Q = 10 m^, 20 mt and ^/s respectively. Here = 172.9 GeV. These results 
show that the renormalization scales for ti production determined by PMC scale setting at 
finite order is insensitive to the choice of very disparate initial scales. 



is slightly different from rrit and its LO cross-section remains almost unchanged. It is noted that 
there is a dip for the NLO scale of the (q'g)-channel, which is caused by the fact that the NLO 
conformal term Aiqq reaches its smallest value and hence the factor B2qq/Aiqq reaches its highest 
negative value when y/s ~ [■\/2exp(5/6)]mt ~ 563 GeV; this results in an exponential suppression 
to the NLO scale. The NLO PMC scale for the (gg)-channel is smaller than rrit by about one order 
of magnitude in low x-region. As a result, its NLO cross-section will be considerably increased; 
i.e. it is a factor of two times larger than its value derived from the conventional scale setting. 
As for the {gg)-ch.aiinel, its NLO PMC scale slightly increases with ^/s, but it is smaller than 
rrit for y/s <C 1 TeV, so that its NLO cross-sections at the Tevatron and LHC are increased by 
15% - 20%. 

There is residual initial renormalization scale dependence because of the unknown-higher-order 
{/3j^^'^}-terms. Because the PMC scales themselves must be a perturbative series of a^, the residual 
initial renormalization scale uncertainty can be greatly suppressed due to the fact that those 
higher-order {/3j^^'^}-terms are absorbed into the PMC scales' higher-order terms. We define a 
ratio -Rn^*" to show how the change of initial renormalization scale affects the PMC scales; i.e. 



,PMC| 

V I 

,,PMC| . 



where /i™*"|^mit=Q stands for the PMC scales determined under the condition of /x™^* = Q, which 
is Ql (LO scale for the non-Coulomb part), Ql* (NLO scale for the non-Coulomb part) or Q2 
(LO scale for the Coulomb part) respectively. In Fig. (14), the ratios for the dominant qq- and gg- 
channels are presented. In order to amplify the differences, we take three disparate scales to draw 
the curves, i.e. Q = 10 nit, 20 nit and a/s respectively. 

As shown in Fig. (14), the LO PMC scale Q2 for the Coulomb-term in both channels are unchanged 
under different choice of Q. Among these choices, Q = usually gives the largest deviation from 
the case oi Q = nit. The residual initial scale dependence for the (5f5f)-channel is small, -R™^ ~ 1, 
only for the LO non-Coulomb PMC scale Ql, it has sizable effect. As an example, for the case of 
Q = y/s, its Ql deviates from that of Q = by ~ 1% at y/s = 7 TeV, and it is raised only up 



56 



Table 4: Dependence on the initial scale of the total tt production cross-sections (in unit: 
pb) at the Tevatron and LHC. Here rrit = 172.9 GeV and the central CTIO as the PDF [164]. 
The number in the parenthesis shows the Monte Carlo uncertainty in the last digit. 





PMC scale setting 


Conventional scale setting 




Q = mt/4: 


Q = mt 


Q = 4:mt 


fir = mt/2 


/ir = rnt 


Hr = 2mt 


Tevatron (1.96 TeV) 
LHC (7 TeV) 
LHC (14 TeV) 


7.620(5) 
171.6(1) 
941.8(8) 


7.626(3) 
171.8(1) 
941.3(5) 


7.623(6) 
171.7(1) 
941.4(8) 


7.742(5) 
168.8(1) 
923.8(7) 


7.489(3) 
164.6(1) 
907.4(4) 


7.199(5) 
157.5(1) 
870.9(6) 



to ~ 7% at ^/s = 14 TeV. In the case of the (gg)-channel, the residual scale dependence of the 
LO/NLO PMC scale for the non-Coulomb part is somewhat larger; i.e. the deviation is about 
12% for the case of Q = ^/s at ^/s = 2 TeV, and the deviation reaches up to ~ 60% at y/s = 14 
TeV. (We expect that this dependence on Q will be greatly reduced at NNNLO.) However in such 
high collision region {^/s > 2 TeV), the total cross-sections are highly suppressed by the parton 
luminosities and their values are almost unchanged by using very different initial scales. 

Total cross-sections with several typical initial renormalization scale /u™'* = Q are presented in 
Table 4. At the NNLO level, it is found that the residual scale uncertainty to the total cross- 
section is less than 10~^ by setting Q = Arrit ot Q = mt/4:. In fact, even by setting Q = 20mt 
and ^/s, such residual scale uncertainty is still less than 10~^ [45]. As a comparison, we also 
present the results for the conventional scale setting in Table 4; by varying the renormalization 
scale within the region of [mt/2, 2 m^], we obtain a large renormalization scale-uncertainty 
at the Tevatron and LHC, which agrees with the previous results derived in the literature, c.f. 
Refs. [146, 147]. This shows that the renormalization scale uncertainty is greatly suppressed and 
essentially eliminated using PMC at the NNLO level. 

We can analyze the combined PDF and as uncertainty by using different CTEQ PDF sets, i.e. 
CTIO [164], which are global fits of experimental data with varying as{mz) G [0.113,0.230]. As 
for the total cross-section after PMC scale setting, we obtain 

C"Tcvatron, 1.96TeV = ^ .626~^qI^I pb (210) 
CTLHC,7TeV = 171.81}^;^ pb (211) 

aLHc,i4Tev = 941.31??;? pb (212) 

where the errors are caused by the PDF+a^ uncertainty. Here a larger PDF+a^ error than that 
of Refs. [147, 149] is due to the choice of PDFs with a wider range of as('^z)- If taking the present 
world average as{mz) — 0.118 ± 0.001 [28], we will obtain a much smaller PDF+a^ error; i.e. 

0"Tcvatron, 1.96TeV = 7.626;i;Q; JgQ pb (213) 
CTLHC,7TeV = 171.8l^;8pb (214) 
t^LHC,14TeV = 941.3ll|^ pb (215) 

The total cross-section ati is sensitive to the top-quark mass, and it is found that the total cross- 
sections decrease with the increment of top-quark mass. After PMC scale setting, by varying 
mt = 172.9 ±1.1 GeV [28], we predict 

CTovatron, 1.96TeV = 7 . 626^^^111 pb (216) 
^LHCTTeV = 171.8l|« pb (217) 
CTLHCMTeV = 941.3l^|4 pb (218) 



57 




Figure 15: Total cross-section au for the top-pair production, with or without PMC scale 
setting, versus top-quark mass. The experimental data are adopted from Refs. [39, 40, 41, 42]. 



where the errors are caused by the top-quark mass uncertainty. In Fig. (15) we present the total 
cross-section function of rrit, where with or without PMC scale setting are shown 

explicitly. After PMC scale setting, the value of au is more closer to the central value of the 
experimental data [39, 40, 41, 42], which shows a better agreement with the experimental data. 

As a summary, after PMC scale setting, we obtain the following points: 

• A larger value for au is obtained, which agrees with the present Tevatron and LHC experimental 
data well. This is achieved because we have resummed the universal and gauge invariant higher- 
order corrections which are associated with the running of the coupling into the LO- and NLO- 
terms by using suitable PMC scales. 

• After PMC scale setting, a more convergent pQCD series expansion is obtained and the resulting 
LO- and NLO- terms are conformally invariant and do not depend on the choice of renormalization 
scheme. The slight change of PMC scales will lead to large effects due to the explicit breaking of 
the conformal invariance [37]. 

• In principle, the PMC scale and the resulting renormalized amplitude is independent of the choice 
of the initial renormalization scale. The residual scale-uncertainty will be greatly suppressed 
when the PMC scales have been set suitably. Then, the usual renormalization scale uncertainty 
Aau/au ~ (-4%) Tevatron and LHC is greatly suppressed or even eliminated by PMC 
scale setting. 



5.2 Top-Quark Pair Backward-Forward Asymmetry 

The top-quark pair forward-backward asymmetry which originates from charge asymmetry physics [165, 
166] has been studied at the Tevatron and LHC. Two options for the asymmetry have been used for 
experimental analysis; i.e. the tt-rest frame asymmetry 

- a{yf>0)-a{yf<0) 
^^^-a(yf>0) + a(yf <0) ^'''^ 

and the pp-laboratory frame asymmetry 

, _ ajyf > 0) - ajyf < 0) 
^^^-a(yf>0) + a(yf <0)' ^''"^ 
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where yf (yf^) is the top quark rapidity in the tt-rest frame (pp-laboratory frame). 

The CDF and DO collaborations have found comparable values in the tt-rest frame: A^^^^ = 
(15.8 ± 7.5)% [167] and A^p^^ = (19.6 ± 6.5)% [168]. The asymmetry in the pp-laboratory frame 
measured by CDF is A^/^^^ = (15.0 ± 5.5)% [167]. The CDF collaboration has also measured the 
dependence of A^p^ with respect to the tt-invariant mass Mtf the asymmetry increases with Mu, and 
A*pg{Mu > 450 GeV) = (47.5 ± 11.4)% [167]. The measured asymmetries are much larger than the 
usual SM estimates for the top quark forward-backward asymmetries. The NLO QCD contributions 
to the asymmetric tt-production yield Ap^ ~ 7% and ~ 5% [169, 170, 171, 172, 173], which are 
about 2(T-deviation from the above measurements. For the case of > 450 GeV, using the MCFM 
program [174], one obtains A*p^{Mit > 450 GeV) ~ 8.8% which is about 3.4o"-deviation from the data. 
A recent reevaluation of the electroweak correction raises the QCD asymmetries by at most 20%: i.e. 
AfsiAfs) ~ 9% (7%) [175, 176] and A%{Mu > 450 GeV) ~ 12.8% [176]. The large discrepancies 
between the SM estimates and the data have aroused interest, because of the possibility for probing 
new physics beyond the SM. However, these comparisons are based on the conventional scale setting 
for choosing the renormalization scale. In the following, we will show how the use of PMC can greatly 
improve our estimates within the SM. 



5.2.1 Basic Formulas 

Numerical results for the top-quark pair production at the Tevatron and LHC have been presented in 
Sec. 5.1. We will compare the total cross-sections derived from PMC scale setting and the conventional 
scale setting. We emphasize two points in addition to the ones listed in the above subsection: 

• At the Tevatron, the top-quark pair production is dominated by the (gg)-channel which provides 
about 85% contribution to the total cross-section. The (gg)-channel due to interference of the one 
gluon and two gluon intermediate states is asymmetric at the NLO level, which leads to sizable 
top-quark forward-backward asymmetry at the Tevatron. In addition, emission of real gluons gives 
an asymmetry. In contrast, one finds that the dominant channel at the LHC is the symmetric 
{gg)-ch.a.nne\, so the top-quark forward-backward asymmetry from other channels will be greatly 
diluted at the LHC; i.e., this asymmetry becomes small which agrees with the CMS and ATLAS 
measurements [177, 178]. 

• More specifically, at the lowest order, the two channels qq — >■ tt and gg ti do not discriminate 
the final top-quark and top-antiquark, so their differential distributions are symmetric for the 
hadronic production. At the NLO level, either the virtual or real gluon emission will cause sizable 
differences between the differential top-quark and top-antiquark production, thus leading to an 
observable top-quark forward-backward asymmetry. At the Tevatron, the asymmetric channels 
are {qq)--, (gq)- and (gq)- channels accordingly. Table 2 shows the total cross-sections of the (gq) 
and (gq) channels are quite small, less than 1% of that of (gg)-channel, so their contributions to 
the asymmetry can be safely neglected. 

Writing the numerator and the denominator of the two asymmetries ApB defined by Eqs. (219,220) 
in powers of a^, we formally obtain 



(221) 







Ni+ as[N2 rr — ] + a' 



y ■ Di Do Do 



where the Dj-terms stand for the total cross-sections at certain a^-order and the A^j-terms stand for the 
asymmetric cross-sections at certain a^-order. The terms up to NLO {Do, Di,Ni) have been calculated, 
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whereas only parts of D2 and N2 are currently known [138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 
148, 149, 150, 151, 152]. 

As shown in Table 2, using conventional scale setting, the relative importance of the denomina- 
tor terms is [a'^Do : a^Di : afD2 ~ 1 : 18% : 12%], and the numerator terms for the asymmetric (qq)- 
channel satisfy [alNi : ~ 1 '■ 50%]. Since at present the NNLO numerator term is not available, 

as a first approximation, we treat these asymmetric terms to have the same relative importance as their 
total cross-sections; i.e. {alNi)qq : {a'^N2)qq ~ (a^Di),^ : {agD2)qq. Thus, the NiDi/Dq term and the 
A^2 term have the same importance. Then, to be consistent, one has to keep only the first term in 
Eq.(221); i.e. dealing with only the so-called LO asymmetry [166, 169, 175, 176]: ApB = ^Q^s- 

On the other hand, after PMC scale setting, we have [a^-Do • ct^-Di : afZ)2 ~ 1 : 41% : 2%] and the 
numerators for the asymmetric (gg)-channel becomes [(y^Ni : a^N2 ~ 1 : 3%]. It shows that, after PMC 
scale setting, the NNLO corrections for both the total cross-sections and the asymmetric part are lowered 
by about one order of magnitude. Therefore, the NNLO-terms and D2 can be safely neglected in 
the calculation, and we can obtain an accurate asymmetry at the NNLO level: 



A - 



Ni-as\ — ] +ots 



Furthermore, it is natural to assume that those higher-order terms Ni and Di with i > 2 after PMC 
scale setting will also give negligible contribution; the above asymmetry can thus be resummed to a 
more convenient form: 

= 2n 1 3n • 222 

Furthermore, as shown in Refs. [166, 175, 176], the electromagnetic and weak contributions A'^o^i provides 
an extra ~ 20% increment for the asymmetry; thus the electromagnetic contribution provides a non- 
negligible fraction of the QCD-based antisymmetric cross-section with the same overall sign. Then, our 
final formula to calculate the asymmetry changes to 

Based on the above considerations, the top-quark forward-backward asymmetry after PMC scale 
setting can be written as 



A 



tt,PMC 



1 



FB tot,PMC (,,VMC\ 



(/^r^; yf > 0) - (/^r^; yf < o)] , (224) 



.pp.PMC _ \ \ (qq) /'..PMC. VP ^ n^l - rr('?'?) (n^^'^- < O^ll r22^1 

^FB — tot,PMc / pMc^ V "''y v"" ) "-^y v"" yJ ' ^ ' 

^ HlH2^tt+X\^^r ) 

where according to Eq.(223): total hadronic cross-section up to NLO; aifj stands for 

the asymmetric cross-section of the (gg)-channel which includes the above mentioned 

and 0{a'^) terms. In the denominator for the total cross-section up to NLO, for each production 

channel, we need to introduce two LO PMC scales which are for the Coulomb part and non-Coulomb 

part accordingly, and one NLO PMC scale for the non-Coulomb part. While in the numerator, we 

only need to know the NLO PMC scale fi^^'^'^^^ for the (gg)-channel, since it is the only asymmetric 

component. 

It is interesting to observe that there is a dip for the NLO scale ^™c,nlo q£ ^^le (gg)-channel when 
y/s ^ [\/2 exp(5 / 6)]mt ~ 563 GeV, which, as shown in Sec. 5.1, is caused by the correlation among the 
PMC coefficients for NLO and NNLO terms. More specifically, it is found that 

^PMCNLO ^ f ^\ ^PMCLO ^ f ^\ g^p f ^ ^ ^init^ ^226) 
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Table 5: Total cross-sections (in unit: pb) for the top-quark pair production at the Tevatron 
with pp-coUision energy -y/i = 1.96 TeV. For conventional scale setting, we set the scale 
^j. = Q. For PMC scale setting, we set the initial scale /i™^* = Q and then apply the 
PMC procedure. Here we take Q = mt = 172.9 GeV and use the MSRT 2004-QED parton 
distributions [179] as the PDF. 





Conventional scale setting 


PMC scale setting 




LO 


NLO 


NNLO 


total 


LO 


NLO 


NNLO 


total 


(gg)-channel 


4.890 


0.963 


0.483 


6.336 


4.748 


1.727 


-0.058 


6.417 


{gg)-Qha,mie\ 


0.526 


0.440 


0.166 


1.132 


0.524 


0.525 


0.160 


1.208 


((7g)-channel 


0.000 


-0.0381 


0.0049 


-0.0332 


0.000 


-0.0381 


0.0049 


-0.0332 


((y'g)-channel 


0.000 


-0.0381 


0.0049 


-0.0332 


0.000 


-0.0381 


0.0049 


-0.0332 


sum 


5.416 


0.985 


0.659 


7.402 


5.272 


2.176 


0.112 


7.559 



where the coefficients are defined through the standard PMC scale setting, cf. Eq.(201). As shown in 
Fig. (12), the value of B2qq is always negative and Aiqq has a minimum value at small y/s. Quantitatively, 
the NLO PMC scale /i™^''^^'-' for the (gg)-channel is considerably smaller than mt in small -y/i-region. 
The NLO cross-section of the (gg)-channel will thus be greatly increased; it is a factor of two times 
larger than its value derived under conventional scale setting, as shown by Table 2. 



5.2.2 Numerical Analysis for the Backward- Forward Asymmetry 

The PMC asymmetries A^g^*" and A^^'^^^ can be compared with the asymmetries calculated using 
conventional scale setting. For definiteness, we apply PMC scale setting to improve HoUik and Pagani's 
results [176], and we obtain 



A 



tt,PMC 
FB 



APP 
^FB 



PMC 



tot.HP 
tot,PMC_ 

tot,HP 
tot,PMC_ 



TyPMCNLO 



HP3 (^conv) 
; (^PMCNLOj 
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a 



«,HP| 
FB I 



0{al) + 



■PMCNLO 



a 



A 



tt,HP| 

FB lOia^a) 



+ A 



ft,HP| 
FB \0 



a 



(/i 



conv j 



^Fif^lci(ai) 



2 (^PMCNLOA 



a 



APP' 



HP I 



FB 10(0^0) 



a^a) + ^^B 




Here stands for the scale set by conventional scale setting and the symbol HP stands for the 

corresponding values of Ref. [176]; i.e. for = mt, it shows [176]: cr^^^^^^tix = 5.621 pb and 



A 



«,HP| 

FB lOia^) 



FB IO(q:3) 



7.32% A^b^|ci(q2q,) - 
:4.85% A^if^oK,) 



1.36% Ap']^^\o{a2) = 

= 0.90% 



0.26% 
= 0.16% 



where 



tf,HP| 

FB 1 0(0!) 



and A^p, 



^FB^^\o{a3) stand for the pure QCD asymmetry at the a^-order under the tt-rest 



A 

frame and the pp lab frame, respectively. 



FB \0{a^a) 
2 
s 



and A^^^^|o(-ct2„) stand for the combined QED and weak with the QCD asymmetry 



A 

at the a^a-order under the tt-rest frame and the pp lab frame, respectively. 



A 



tt,HP| 

FB IO(a2) 



2\ and 



HP I 



FB 



stand for the pure electroweak asymmetry at the a^-order under the 



tt-rest frame and the pp lab frame, respectively. 
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Figure 16: Dominant cut diagrams for the n/-terms at the a^-order of the (gg)-channel, 
which are responsible for the smaller effective NLO PMC scale /I™*-'''^'"^, where the solid 
circles stand for the light-quark loops. 

0.35 r 




Figure 17: Comparison of the PMC prediction with the CDF data [167] for the tt-pair 
forward-backward asymmetry for the whole phase-space. The Hollik and Pagani's results 
(HP) [176] using conventional scale setting are presented for a comparison. The result for 
DO data [168] shows a similar behavior. 



Total cross-sections for the top-quark pair production at the Tevatron with pp-coUision energy a/s = 
1.96 TeV and with the same parameters of Ref. [176] are given in Table 5. In the formulas (227,228), 
we have defined an effective coupling for the asymmetric part, which is the weighted 

average of the QCD coupling for the (gg)-channel; i.e. in using the effective coupling (jJr^'"'^^'^^ , 
one obtains the same (gg)-channel NLO cross-section as that of as(/i™*"''^'"^)^. 

It is noted that the NLO-level asymmetric part for (gg)-channel only involves the NLO PMC scale for 
the non-Coulomb part, so the effective coupling Og (^/I™^''^^*^) can be unambiguously determined. We 
obtain a smaller effective NLO PMC scale ^™c,effective ^ exp(— 9/10)mt ~ 70 GeV, which corresponds 
to as (/I™^'^^^) = 0.1228. It is larger than af ^ (mt) ^ 0.098 [175, 176]. This effective NLO PMC 
scale is dominated by the non-Coulomb rij-terms at the a^-order, which are shown in Fig.(16). In these 
diagrams, the momentum flow in the virtual gluons possess a large range of virtualities. This effect for 
NLO PMC scale ^™c,effectivc regarded as a weighted average of these different momentum flows 

in the gluons, which can be softer than the nominal scale, rrit. Finally, we obtain 

^%^^~12.7%; Afij™^ ^ 8.39% . (229) 



^In principle, one could divide the cross-sections into symmetric and asymmetric components and find PMC scales 
for each of them. For this purpose, one needs to identify the Tiy-terms or the n^-terms for both the symmetric and 
asymmetric parts at the NNLO level separately. 
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Thus after PMC scale setting, the top-quark asymmetry under the conventional scale setting is in- 
creased by ~ 42% for both the tt-rest frame and the pp-laboratory frame. This large improvement is 
explicitly shown in Fig. (17), where Hollik and Pagani's results which are derived under conventional 
scale setting [176] are presented for comparison. 

Another possible effect from QCD can be the lensing effect of the final state interactions of the t and 
t with the beam spectators. The same diagrams causes Sivers single-spin asymmetry and diffractive 
deep inelastic scattering^. 





A^g (M^>450GeV) 






±2a 






HP 

PMC 

























Figure 18: The PMC prediction of Ap^{Mtt > 450 GeV) and the corresponding CDF 
data [167] for the tt-pair forward-backward asymmetry for Mu > 450 GeV. The Hollik and 
Pagani's results (HP) [176] using conventional scale setting are presented for a comparison. 



The CDF collaboration has found that when the tt-invariant mass, Mu > 450 GeV, the top- 
quark forward-backward asymmetry A^p^lMu > 450 GeV) is about 3.4 standard deviations above 
the SM asymmetry prediction under the conventional scale setting [174]. After applying PMC scale 
setting, we have a°l'™_^^^xiMtt > 450 GeV) = 2.406 pb and a, (71™^'^^°) = 0.1460 with 71^^^'^^° ~ 
exp(— 19/10)mt ~ 26 GeV. Then, we obtain 

A%^^^{Mu > 450 GeV) ~ 35.0% , (230) 

which is increased by about 1.7 times of the previous one A'Pb {Mtt > 450 GeV) = 12.8% [176]. Our 
present prediction is only about la-deviation from the CDF data, which is shown in Fig. (18). This 
shows that, after PMC scale setting, the discrepancies between the SM estimate and the present CDF 
and DO data are greatly reduced. 



6 Summary 

Because of the RG invariance (38,39), the predictions for a physical observable must be independent 
of the renormalization scheme and the initial scale. The results cannot depend on which scheme the 
theorist chooses; e.g. MS'-scheme, MOM-scheme, etc. Note that the conventional M5'-scheme is 
somewhat artificial. One can introduce a more general MS'-like renormalization scheme, T^-^-scheme, 
by further absorbing an arbitrary constant 6 into 1/e pole, i.e. ^ + ln(47r) — 7^; — 5. Physical results 
cannot depend on the choice of S. 

^We thanks Benjamin von Harling and Yue Zhao for conversions on this possibihty. 
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At a fixed-order the dependence on the renormahzation scheme and initial scale choice leads to large 
uncertainties for perturbative QCD predictions. The problem is compounded in multi-scale processes. 
The conventional scale setting procedure assigns an arbitrary range and an arbitrary systematic error 
to fixed-order pQCD predictions. As we have discussed in this review, this ad hoc assignment of the 
range and associated systematic error is unnecessary and can be eliminated by a proper scale setting 
as the PMC. 

The extended RG equations, which includes the dependence on the scheme parameters, provide a 
convenient way for estimating both the scheme and scale dependence of the perturbative predictions 
for a physical process. It provides a way for the running coupling to run reliably either in scale or in 
scheme. With the help of the extended RG equations, we have presented a general demonstration for 
the RG invariance. Furthermore, this formalism provides a platform for a reliable error analysis, and it 
also provides a precise definition for the QCD asymptotic scale under any renormahzation 7^-scheme, 
^QCD^^ which is defined as the pole of the strong coupling in the 't Hooft scheme associated with 
7?.- scheme. 

Several scale setting methods have been proposed in the literature: FAC, PMS, BLM and PMC. 
The FAC (Fastest Apparent Convergence) use the scale to contract the prediction to one term. The 
PMS (Principle of Minimum Sensitivity) chooses the scale at the point of minimum variation. The 
BLM (Brodsky-Lepage-Mackenzie) and PMC (Principle of Maximum Conformality) procedures shift 
all {/3j}-terms into the argument of the running coupling. Based on the extended RG equation, we 
have discussed the self-consistency conditions for a scale setting method, which include the existence 
and uniqueness of the renormahzation scale, reflexivity, symmetry, and transitivity. These properties 
are natural requirements of RG invariance. We have shown that the FAC and BLM/PMC satisfy 
these requirements, whereas the PMS does not. The PMS is designed to be renormalization-scheme 
independent; however it violates the symmetry and transitivity properties of the renormahzation group, 
and does not reproduce the Cell Mann-Low scale for QED observables. In addition, the application 
of PMS to jet production from the e+e~ — > qqg gives unphysical results; i.e. its PMS scale rises 
without bound for small jet energy, since it sums physics into the running coupling not associated 
with renormahzation. This implies the necessity of further careful studies of the theoretical principles 
underlying PMS. 

Among these scale setting methods, the advantages of PMC are clear. In PMC, the same procedure is 
valid for both space-like and time-like arguments; in particular, this leads to a well-behaved perturbative 
expansion, since all the large {/3i}-dependent terms on the time-like side involving vr^-terms are fully 
absorbed into the running coupling. Through the PMC - BLM correspondence, the PMC and the BLM 
are equivalent to each other. Thus the features of BLM scale setting are also adaptable to PMC. 

For convenience, we summarize the dominant features of PMC in the following: 

• It keeps the information of the higher order corrections but in a more convergent perturbative 
series. After PMC scale setting, the divergent "renormalon" series with n!-growth disappear, so 
that a more convergent perturbative series is obtained. 

• Its estimation is renormalization-scheme independent, because after PMC scale setting, 

— the resulting expressions are conformally invariant and thus do not depend on the choice of 
renormahzation scheme; 

— one obtains the proper scale-displacements among the PMC scales derived for different 
schemes or conventions; 

— one also obtains GSRs connecting observables and schemes. 

These features become clear using the T^^-scheme. The 5-terms always accompany nonconformal 
/9-dependent terms, and thus the elimination of 5-terms by shifting the scale of the running 
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coupling is equivalent to the elimination of {/3i}-terms. The PMC estimate can therefore also be 
achieved through a proper treatment of (5-terms. This new way for PMC scale setting can be 
readily programmed for automatically setting the PMC scales to all orders [48, 129]. 

• The PMC provides a systematic way to set the optimized renormalization scale for a fixed-order 
calculation. In principle, the PMC needs an initial value to initialize renormalization scale. It is 
found that the estimates of PMC are to high accuracy independent of the initial renormalization 
scale; even the PMC scales themselves are in effect independent of the initial renormalization scale 
and are 'physical' at any fixed order. This is because the PMC scale itself is a perturbative series 
and those unknown higher-order {/3j}-terms will be absorbed into the higher-order terms of the 
PMC scale, which is strongly power suppressed. 

• By applying PMC to the top-quark pair hadroproduction up to NNLO level, it has been found 
that the PMC scales and the resulting finite-order total cross-sections are both to high accuracy 
independent of the choice of an initial scale [37, 45, 127]. After PMC scale setting, the top 
quark forward-backward asymmetries at the Tevatron are : A^p^ ^ 12.5%, A^p^ ^ 8.28% and 
A*p^{Mtt > 450 GeV) ~ 35.0% [127]. These predictions deviate approximately 1 a from the CDF 
and DO measurements. The large discrepancy of the top quark forward-backward asymmetry 
between the standard model estimate and the data is thus greatly reduced. 

• A PMC analysis for the charmonium production processes, + e~ J/tp{tlj') + Xcj with (J = 
0, 1, 2), has been given in Ref. [180]; it shows that the scale uncertainty for both the polarized and 
the unpolarized cross sections are greatly suppressed even at the NLO level. 

• The PMC is adoptable for the QED case. In the Abelian limit Nq — ?■ at fixed a = CpUg with 
Cp = (N^ — i)/'2Nc [29, 36], the PMC also agrees with the standard Cell Mann-Low procedure 
for setting the renormalization scale in QED. Any method used in perturbative QCD must be 
applicable to QED in the Nc — ?■ limit, since they share the same Yang-Mills Lagrangian. 

The elimination of the renormalization scale ambiguity and the scheme dependence using the PMC 
procedure will not only increase the precision of QCD tests, but it will also increase the sensitivity of 
collider experiments to new physics beyond the SM. The PMC procedure can be advantageously applied 
to the entire range of perturbatively-calculable QCD and Standard Model processes, eliminating an 
unnecessary systematic error. 
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